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Abstract 



The usefulness of the Relativistic Schrodinger Theory (RST) is studied in the 
field of atomic physics. As a concrete demonstration, the positronium ground- 
state is considered in great detail; especially the groundstate energy E is worked 
out in the non-relativistic approximation and under neglection of the magnetic 
interactions between the positron and the electron. The corresponding RST 
prediction (Eq ~ — 6,48[el / ]) misses the analogous conventional Schrodinger re- 
sult (E ~ —6,80 [eV]) but is closer to the latter than the corresponding Hartree 
approximation (—2,65 [eV]). The missing binding energy of 6, 80— 6, 48 = 0,32[eV] 
can be attributed to the approximative use of an SO (3) symmetric interaction po- 
tential which in RST, however, is actually only SO (2) invariant against rotations 
around the z-axis. It is expected that, with the correct use of an anisotropic inter- 
action potential due to the SO (2) symmetry, the RST predictions will come even 
closer to the conventional Schrodinger result, where however the mathematical 
structure of RST relies on exotic (i.e. double-valued) wave functions and on the 
corresponding unconventional interaction potentials (e.g. Struve-Neumann poten- 
tial). 
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I. INTRODUCTION AND SURVEY OF RESULTS 



The general philosophy, being well accepted by most of the physicists, says 
that any new theory must be in agreement with the successes of the existing the- 
ories of long standing. There are many instances in the development of physics 
which can serve as a demonstration of this dogma: general relativity vs. Newto- 
nian gravitation theory , special relativity vs. Newtonian mechanics |3|, |4|, 
Maxwellian electrodynamics vs. Faraday's theory of magnetism and Gaufi' law of 
electrostatics [5|, or Heisenberg's matrix mechanics and Schrodinger's wave me- 



chanics vs. Bohr's and Sommerfeld's old quantum theory [(3] etc. The present 
paper is intended to present a further example of this hierarchical "inclusion" of 
physica 
(RST) 

relativistic Schrodinger quantum mechanics for many-particles systems, must be 
required to reproduce the well-known and generally accepted results of the con- 
ventional quantum mechanics! Clearly, such a great goal can hardly be achieved 
in full generality by one single paper; but for a first demonstration of this one 
would be satisfied with a special but nevertheless typical situation so that there 
may arise some confidence into the validity of the inclusion dogma for the general 
case of RST. 

Our choice of such a demonstration refers to the binding energy (Eq) of the 
positronium groundstate in the non-relativistic approximation and with neglection 
of the magnetic (i.e. spin-spin) interactions between the electron and the positron. 
In this lowest order of approximation, the conventional Schrodinger theory predicts 
the binding energy E (= -6,80 eV) by means of a very simple argument [10]: The 
two-particle wave function ^(fi,f 2 ) is factorized into the product ip{R) ■ ip(r) 
when written down in terms of the center-of-mass coordinate R (= | • (ri + r*2)) 
and relative coordinate r (= (fx — F2)); and then the conventional two-particle 
Schrodinger equation (see equation (16 . 1 5f) below) yields the ordinary Coulomb 



theories; namely the recently established Relativistic Schrodinger Theory 
as a non-standard relativistic generalization of the ordinary non- 



force problem for the internal motion written in terms of the relative coordinate r. 
However the solution of this standard problem in atomic physics is well-known and 
the conventional groundstate energy for positronium is easily found as 

£o | conv . = -~ = ~\<4Mc 2 - -6,80 [eV] , (1.1) 

fc2 2 

where cl-q is the Bohr radius (= -S-?) and «g(= jb) is the fine structure constant. 
Indeed this result is nothing else than the hydrogen groundstate energy due to 
a fixed nucleus, with merely the electron mass M being replaced by the reduced 
mass 4f on account of the comoving positron. Clearly there exists an extender 



Uterature for i— g these se.eielasacal predictions for both the hydrogen PJ| 
and positronium [121] semiclassical groundstate energy Eq, namely by taking into 
account various kinds of QED corrections. However the point here is that one 
would like to see first the coincidence of the results in the lowest-order of approxi- 
mation before one compares the higher-order predictions of the various theoretical 
approaches. This is the reason why the intention of the present paper must consist 
in elaborating (within the RST formalism) the positronium groundstate energy Eq 
in the non-relativistic approximation and with neglection of the spin! According 
to the inclusion dogma, this RST result should then agree with the conventional 
prediction (11. ip 

However the calculation of the non-relativistic Eq is within the RST formalism 
not so simple as in the conventional approach since both kinds of approximations 



do require some su 
ceding papers 



e considerations: First, it has been asserted in some pre- 
14j that the non-relativistic limit of RST coincides with the 
well-known Hartree-Fock approach. However the problem is here that the latter 
approach is an approximation method whose predictions do not always come close 
enough to the corresponding predictions of the standard Schrodinger quantum 
mechanics. A drastic example for this deficiency refers just to the positronium 
groundstate Eq which is missed by the Hartree approximation by roughly 50%! 
(see Sect. VI. B). Thus if, for any physical situation, the non-relativistic limit of 



RST would agree with the Hartree-Fock approach, RST would be afflicted with 
the same deficiencies as HF; and it would therefore probably not be worth while to 
further elaborate the new formalism. However the subsequent RST treatment of 
positronium demonstrates the occurrence of "exotic" states and the correspond- 
ing interaction potentials which shift the corresponding RST predictions closer 
towards the conventional Schrodinger results (and thus also closer to the observa- 
tional data) than it is possible for the original Hartree-Fock method (Sect. VI). 

A second point of concern refers to the spherical symmetry of the (internal) 
groundstate of the conventional approach. This spherical symmetry is of course 
due to the corresponding symmetry of the (internal) Coulomb potential. However, 
if one approaches the non-relativistic limit of RST from the side of the fully rela- 
tivistic formalism, it is not a matter of course that one will encounter a spherically 
symmetric interaction potential for the two particles. The reason for this is that 
the relativistic approach takes account of the particle spin (~-> Dirac equation); and 
usually there emerges some preferred direction in three-space when angular mo- 
mentum (spin or orbital) comes into play in connection with the bound states 
angular momentum conservation). Especially for the present situation with RST, 
the preferred space direction leaves its imprint on the electric interaction potential 
in such a way that its spin-induced anisotropy survives the non-relativistic limit. 
Therefore, even in this limit, we have to deal with an anisotropic potential (see the 
Struve-Neumann potential fl4.17l) - fl4.18l) ); but we will neglect this anisotropy and 
will take into account only the spherically symmetric part in order that our non- 
relativistic treatment displays the same symmetry as the conventional theory (i.e. 
the spherical symmetry of the internal Coulomb potential). Of course it should be 
clear that one has to pay some prize for such a forceful approximation procedure, 
and this means concretely that our non-relativistic RST prediction for the positro- 
nium groundstate energy E will adopt a certain inaccuracy (see the discussion 
of this point in the Appendix). Actually, in place of the expected 6,80 [eV] of 



5 



equation (II. ip . one ultimately arrives at 6,48 [eV] in the spherically symmetric 
approximation, so that the elaboration of the proper RST prediction will require a 
more efficient approximation technique for anisotropic potentials (not to be treated 
in the present paper). 

The results, being described so far, are worked out by the following arrangement: 

A. Relativistic Schrodinger Theory 

In Sect. II, the general RST dynamics for two-particle systems is briefly 
sketched in order that the paper be sufficiently self-contained (for more details 



le paper b 



see the preceding papers |7|]-|9|], Il3]-[l5]). Here the key points refer to the cou- 
pled matter and gauge field equations (i.e. Relativistic Schrodinger Equation ( 12.15P 
and non-Abelian Maxwell Equation (12.381) ). and also to the specific gauge struc- 
ture of the theory which admits to treat the electromagnetic and exchange forces 
on the same footing. Furthermore the general RST formalism provides one with 
the energy- momentum densities T^ u (e.g. for matter, see (I2.2ip below), the time 
components (Too) of which serve later on to define the energy Et (I5.2aj) of the 
bound stationary field configurations. 

B. Non- Standard Eigenvalue Problem 

Sect. Ill presents the mass eigenvalue equations for the stationary two-fermion 
systems. The crucial point here is that, at least for the ground state, both particles 
must occupy physically equivalent states which can differ at most with respect 
to the relative orientation of their magnetic fields H a (f),a = 1,2. This simple 
logical requirement leads to two classes of possible field configurations: ortho- 
positronium (13. 7aj) and para-positronium (13. 7b[) . Naturally the two wave func- 
tions ip a (r) (a = 1, 2) due to both classes turn out to be very similar but are not 
identical as long as the magnetic (i.e. spin-spin) interactions of both particles are 
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taken into account, see the eigenvalue equations for ortho-positronium (I3.29al) - 
(13.29bl) and for para-positronium (I3.32ap - (13.32bj) . If, however, the magnetic in- 



teractions are neglected, the spin-up and spin-down components of the wave 
functions ip a {r) become decoupled and this leads to a simpler eigenvalue sys- 
tem, (I3.57al) -( l3.57dl) . which has the advantage to admit spherically symmetric so- 
lutions. Their non-relativistic limit can easily be obtained, cf. (I3.66al) or (14 .7p . 
and does serve later on as a trial function for minimalizing the groundstate energy 
functional (Sect. VI). 

The most important point, however, refers here to the fact that the solu- 
tions <p±(r,9,(fi) of the mass eigenvalue system are double- valued [9|], cf. (13. 55[) : 
and their non-relativistic limit R (I3.66aj) does not obey an ordinary one-particle 
Schrodinger equation (as is usually the case with the Hartree-Fock approxima- 
tions) but it obeys a modified form thereof, see equation (13.681) . This is one of 
the origins of the potentiality of RST to predict the energy of the bound systems 
in better agreement with the experimental data (Sect. VI) than is possible for the 
Hartree-Fock approach! 

C. Unconventional Gauge Potentials 

Sect. IV reveals the second origin of the superiority of RST over the Hartree- 
Fock approach: this refers to the gauge potential which mediates the interaction 
between the two particles. Since the magnetic interactions are neglected, one is 
concerned here exclusively with the time component Ao(r) (of A^) and this is 
linked to the charge density jo(r) of the particles by a Poisson equation, cf. (I3.18P 
below, which itself is to be deduced from the Maxwell equation (12.381) . However the 
point here is that the RST ground state of positronium is of a rather exotic nature, 
namely in the sense that the bosonic character of the bound two-fermion system 
is transferred to anyone of its constituents: they individually own an integral 
quantum number of total angular momentum (J z = 0), see equation (13.411) . 
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A consequence of this exotic behavior of the fermions refers to the electrostatic 
potential A Q (r) (14.81) being generated by the corresponding electric charge distribu- 
tion jo{r). On principle, such an exotic potential A${r) must always be anisotropic 
(even for the groundstate), but in the Appendix it is demonstrated that for a first 
rough estimate of the positronium groundstate energy Eq it is sufficient to use 
only the isotropic part ^A (r) of the electrostatic interaction potential. This 
spherically symmetric part ^Ao(r) is defined by equation (14.111) ; and when the 
simple exponential function R(r) (14. 7p is taken as the (non-relativistic) trial wave 
function, the associated potential is the Struve- Neumann potential ^Aq(t) ( 14.171) - 
( I4.18P being exhibited in Fig. 1. In contrast to their Hartree-Fock counterparts, 
these exotic RST potentials are still singular at the origin (r = 0), but less irregular 
as the Coulomb potential (~ r -1 ), so that their energy content is finite. 

D. Groundstate Energy 

Concerning the physical relevance of any new theory of quantum matter, per- 
haps the most important object is here the energy functional (E^) whose value 
upon the solutions of the eigenvalue problem yields the energy being carried by 
the bound stationary field configurations. Clearly, the reason is that the set of 
these energy values constitutes the theoretical level system which may then be 
opposed to the observational data in order to test the usefulness of the theory. 
Therefore setting up the right energy functional is one of the crucial points of 
any particle theory. In this sense, Sect. V presents a thorough inspection of the 
RST energy functional E?, especially concerning its non-relativistic approxima- 
tion, because the non-relativistic groundstate energy of positronium is the main 
concern of the present paper. According to the general expectation, the desired 
energy functional (Et) should be equipped with the property that its value upon 
the groundstate solution of the eigenvalue problem yields the minimally possible 
energy value (E , say) in comparison to all other field configurations. In ordi- 
nary, non-relativistic quantum mechanics such a principle of minimal energy is 
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well-known, i.e. the Ritz variational principle 1 171 ]. whose associated Euler- 
Lagrange equation is just the conventional Schrodinger equation. 

Thus the task is now to find first the right relativistic generalization of the 
Ritz principle within the RST framework; and afterwards one may look for its 
non-relativistic approximation in order to calculate with its help the positronium 
groundstate energy. But here the interesting point is now that the non-relativistic 
approximation of the RST principle of minimal energy does not lead us back to its 
philosophical origin (i.e. the Ritz variational principle); rather, a new form of varia- 
tional principle is found which yields its associated groundstate energy closer to the 
experimental data than is obtained by the well-known Hartree (or Hartree-Fock) 
approach. This RST result is obtained by being satisfied with the spherically sym- 
metric approximation of the groundstate wave functions (Appendix). It remains 
to be tested whether or not the present RST principle of minimal energy would 
yield even exact coincidence with the conventional Schrodinger energy eigenvalue 
if the anisotropy of the RST groundstate is taken into account. 

The RST principle of minimal energy is constructed by the following steps: 
First, the coupled RST system of eigenvalue and gauge field equations is identified 
as the set of Euler-Lagrange equations being due to the well-known Hamilton- 
Lagrange action principle, see equations (I5.8j) - (l5.1ip . For this variational proce- 
dure, the imposition of some constraint is not necessary, neither for the wave func- 
tion ^ nor for the gauge field A^. However the stationarity of the action PVrst with 
respect to the scaling transformations of the electric gauge potential A yields a 
certain relationship between the latter potential and the matter field (i.e. the Pois- 
son identity f )5.20p ) which plays an important part for the subsequent construction 
of the RST principle of minimal energy! Next, the mass functional M^c 2 f)5.2ip 
is constructed as an intermediate step, where this functional is to be considered 
as the immediate relativistic generalization of the conventional Ritz variational 
principle. Here, the first constraint comes into play (i.e. the normalization of the 
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wave functions, see equation (15.221) ). 

Although the mass functional correctly reproduces the mass eigenvalue system 
in form of the corresponding Euler-Lagrange equations, it nevertheless cannot be 
identified with the wanted energy functional (E^) because the gauge field equa- 
tions cannot be deduced from it by means of the usual variational procedure. The 
solution of the problem is finally attained by reconsidering the total field energy E^ 
as the spatial integral of the total energy- momentum density ( T -*7oo(r), see equa- 
tion (I5.2al) . Indeed if the stationary field configurations of this functional are 
determined, under both constraints of wave function normalization and Poisson 
identity, then the coupled system of mass eigenvalue plus gauge field equations is 
properly recovered for the stationary bound systems. Thus, incorporating the two 
constraints by the method of Lagrangian multipliers, the final form of the desired 
RST energy functional E T is given by equation (15.431) . 

E. Numerical Test 

Clearly, the usefulness of any theoretical construction can be evaluated only 
by opposing its predictions to the observational data. Therefore the value of the 
obtained energy functional E^ (I5.43j) upon the chosen trial function R(r) ( \4.7h for 
the positronium groundstate is calculated in Sect. VI and the result is compared 
to the analogous procedure for the Hartree-Schrodinger approximation method 
(since the spin is neglected here, the Hartree-Fock approach reduces to the much 
more simpler Hartree approximation method). The difference between RST and 
the Hartree-Schrodinger approach lies mainly in the interaction potential of the 
electron and positron (see equations (I4.17I) -( 14TT8|) and (\6.27h together with fig. 2) 
where the competition ends in favour of RST: the RST predictions for the positro- 
nium groundstate binding energy (— E ) amount to 6,48 [eV] and therefore are 
closer to the experimental value (6,80 [eV]) than the predictions of the simple 
Hartree-Schrodinger approach (2,65 [eV]), see fig. 2. Moreover it is demonstrated 
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in the Appendix that the residual RST inaccuracy of (6,80 - 6,48) = 0,32 [eV] is 
due to the neglected anisotropy of the Struve-Neumann potential ^A (r), so that 
one may legitimately expect a higher degree of coincidence between the RST pre- 
dictions and the experimental data as soon as the RST approximation techniques 
become improved. 



II. TWO-FERMION SYSTEMS IN RST 



RST, in general, provides a re 
of positive and negative charges 



ativistic description for an arbitrary assemblage 
; but for the present treatment of positronium 



one is concerned with one positively and one negatively charged Dirac particle. 



A. Conservation Laws 



Therefore the total velocity operator F M for such a two-fermion system must 
appear as the direct sum of two ordinary Dirac matrices 7^: 

r M = (-7 M )©7M, (2.1) 

where the minus sign for the first particle (i.e. the positron) is merely due to 
convention (the preceding papers were concerned with N-electron systems for which 
the total velocity operator has been adopted as T M = 7^ © 7^ © . . . ). With the 
existence of such a velocity operator (12 .ip it becomes possible to associate a 
conserved current j M with any two-particle wave function ty(x) which itself is the 
direct sum (Whitney sum) of the single-particle wave functions ip a (a = 1, 2) 

* = Vi © ^2 , (2.2) 

namely by means of the following prescription 

3fl = ^T^. (2.3) 
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Thus the desired conservation law may be written down in local form as 

V% = , (2.4) 
and the associated total charge can be defined through 

z= [ J»dS» (2.5) 



(5) 

where the space-like three-dimensional hypersurface (S) is arbitrary just on ac- 
count of the local law ( 12.41) . Clearly, for the considered positronium, the total 
charge z (12.51) is zero because both particles (i.e. electron and positron) carry 
opposite charges. 

Besides the charge conservation law (12.41) . there are further conservation laws 
which have to be obeyed by a closed two-particle system, e.g. the energy- 
momentum conservation law 

V" (T) 7^ = (2.6) 

where ^ T 'T IM , is the total energy-momentum density of the two-particle system. 
Clearly when this system is non-closed in the sense that it is acted upon by an 
external field ^'F^ generated by the external four-current < - ex ^ M , then the source of 
( T ^T r ,, (\2.^ becomes non-zero and actually is determined by the well-known Lorentz 
force ( xe 'f v , i.e. 

V^ T % u = -^f v (2.7a) 
{xe) fv = -hcF,„-^f ■ (2.7b) 

According to Newton's actio = reactio, this is the reactive force to that force by 
which the total field of the two-particle system pulls at the external current 
(ea % (see the discussion of this point in the preceding papers 

0-0) 

. However, 

the presently considered positronium is assumed to be closed with no external 
field being present (i.e. ^ ex F^ v = 0; ( - ex ^ M = 0), so that the energy- momentum 
conservation (12.61) is strictly valid. 

12 



The problem is now to define the total energy-momentum density ^ T ^T^ U of the 
two-particle system in such a way that the required conservation law ( 12. 6ft . and 
also the charge conservation (12.41) . does automatically hold as an implication of 
the RST dynamics to be readily set up. The solution to this problem is obtained 
by first building up the total density ^^T^ from a matter part ^^T^ and a gauge 
field part ^ G ^T^, i.e. 

{T %„ = {D %„ + , (2.8) 

and then one introduces an energy-momentum operator 7^, for the matter part 
which serves to define the matter density ^ D ^T^ through 

(°% u = (2.9) 

quite analogously to the procedure with the total current j M (12.31) . The problem 
has thus been deferred from the density ^^T^ to the operator 7^„ and the gauge 
field density ^ G ^T Ml/ (see below); but it can be demonstrated that the right choice 
for the matter operator is the following 

% v = -{r^Hv + H U T^ + F v Hp + H^} . (2.10) 

Here a new object is introduced, i.e. the Hamiltonian H^, which will require its 
own field equations. 



B. Hamiltonian Dynamics 

The RST dynamics must be set up now in such a way that the above mentioned 
conservation laws are automatically implied. Evidently, the RST kinematics refers 
to four fields: total velocity operator IT M (12. ip . wave function \I/ (12.21) . gauge field 
(see, e.g., (12.7b[) ) . and Hamiltonian Ti^ (12.101) . Thus, anyone of these objects 
requires the specification of a field equation, and the totality of these field equations 
will then constitute the RST dynamics. 
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Concerning the Hamiltonian H.^, its dynamical equations consist of the inte- 
grability condition 

V^H U - VJHp + ^ \H» H„] = i UcT^ (2.11) 

and of the conservation equation 

V^-^H^ = -ihe{(^y + , (2.12) 

where M. is the mass operator and are the Spin(l, 3) generators constructed 
by means of the velocity operator F M as 

E^ = ^[r M ,rJ . (2.13) 

It can be shown that the Hamiltonian dynamics ( I2.11l) -( l2.12p guarantees the va- 
lidity of both the bundle identities such as, e.g., for the wave function ^ 



T> V — V V 



* = ^> • * (2.14) 



and also the validity of the desired conservation laws (12.41) and ( 12. 6ft . But clearly, 
for the latter purpose one has to specify also an appropriate field equation for the 
matter field 

C. Matter Dynamics 

The field equation for the wave function \l/ is adopted as the Relativistic 
Schrddinger Equation (RSE) 

ihcV^ = H^. (2.15) 

Differentiating this equation once more and applying also the conservation equa- 
tion (l2~T2l recasts the RSE (EHo) to the Klein Gordon Equation (KGE) 

Vn>^ + \ -r-) * = ~^ V F^ • (2.16) 
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However, there is another nice possibility in order to eliminate the Hamiltonian 
Ti^ from the RSE (j2.15p . Namely, one can show that the differential form (I2.12p 
of the conservation equation is equivalent to the following algebraic formulation: 

Mc 2 = WT^ = T^n„ . (2.17) 

From here it may be seen that both the mass operator Ai and the velocity op- 
erator F^ are assumed to be Hermitian (i.e. M. = Ai, F M = T^), but not the 
Hamiltonian Ti^ (7^ 7Y M ). But now multiply through the RSE (12. 15ft from the left 
by the velocity operator F^ and use the algebraic conservation equation (I2.17P in 



order to find the Dirac Equation (DE) for the two-particle wave funcion 

ihT fl V^ = Mc^. (2.18) 

However once the DE is established, it is very instructive to verify the de- 
sired charge conservation (I2.4p just with its help (the verification of the energy- 
momentum conservation (12. 6ft is a little bit more complicated, see e.g. ref.|l5|). 
Indeed, carry through the differentiation process in equation (12.41) and find with 
the help of just the DE (12.181) that the charge conservation actually holds, provided 
the velocity operator F^ is covariantly constant: 

2VT„ = . (2.19) 

In a similar way, one verifies the energy-momentum conservation (12. 6p under the 
additional condition of covariant constancy of the mass operator Ai 

XyW = . (2.20) 

These constancy conditions do attribute the status of absolute objects to F M and 
Ai which however is a less dramatic insight because the Dirac matrices 7 M and 
the particle rest masses M are always taken as absolute objects in most of the 
modern theories of elementary particles. Finally observe also that the DE ( 12 . 181) 
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can be used for the elimination of the Hamiltonian Ji^ from the matter density 
{D} T^ fl23D- fl2~TUD which then adopts the following shape 

{D % V = — {*r M (X> v *) - {V^)T^ + W V (V^) - (Vfi)F„*} . (2.21) 

This form will subsequently be used for the calculation of the positronium ground 
state energy. 

However, what is a really interesting point with those constancy conditions 
is the fact that they reduce the structure group of the theory for non-identical 
particles, whereas no reduction occurs for identical particles! In physical terms, 
this reduction implies the elimination of the exchange forces between the non- 
identical particles, whereas the identical particles undergo both the conventional 
gauge forces (here: electromagnetic interactions) and the exchange forces! This 
effect will readily become elucidated by inspecting now the gauge field equations. 



D. Electromagnetic and Exchange Forces 

An important question in connection with the positronium system refers to 
the phenomenon of exchange forces between the electron and the positron. In 
the conventional theory [lOj, the Schrodinger-Hamiltonian H{f\,r2) is symmetric 
under exchange of the one-particle variables (fi,pi) <-> (r*2,p2) 

and therefore the energy eigenfuntions ^(fi,/^) should be either symmetric or 
antisymmetric under particle exchange (fi r 2 ). Properly speaking, such a 
symmetry of the energy eigenfunctions signals the presence of exchange forces (in 
addition to the usual electromagnetic interactions) with the corresponding shift of 
the energy levels due to the exchange energy. 

On the other hand, the conventional non-relativistic treatment of positronium 
makes use of the transformation from the single-particle coordinates ri,r 2 to the 
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center of mass coordinate R and relative coordinate r, and then there arises a 
one-particle eigenvalue problem for the bound states being written in terms of 
the relative coordinate f alone. The corresponding positronium spectrum emerges 
thus as the usual hydrogen level system (Coulomb force problem) with merely the 
electron mass M being substituted by the reduced mass M/2. It should be imme- 
diately evident that for such an effective one-particle problem there cannot emerge 
exchange effects because the latter do always refer to many-particle systems! For 
the subsequent RST treatment of the positronium spectrum, the exchange forces 
between the electron and positron are also suppressed but on account of quite dif- 
ferent reasons. Therefore it is very instructive to first inspect this RST mechanism 
of suppression of the exchange interactions. 

Observe also that the conventional Hamiltonian H (12.221) works with the exact 
Coulomb force between both particles, with no shielding effects included. In con- 
trast to this, the subsequent RST treatment will lead to some weakening (but not 
complete regularization) of the Coulomb singularity. 

The coupling of the RST matter field $ to the gauge field (bundle connec- 
tion), with associated curvature J 7 ^ 

•F = V fx*Au + [<A.fi, A u ] j (2.23) 

occurs as usual via the principle of minimal coupling; i.e. the gauge- covariant 
derivative of $ in the RSE (f2TT5l) or DE (f2~T8l) is written as 

= + A^ . (2.24) 

Since the two-particle wave function \1/ is here the (Whitney) sum of the two single- 
particle wave functions ip a (a = 1,2), the abstract derivative (12.241) may be recast 
into component form and then looks as follows [8j]: 

= ' i A \^i ~ i (2.25a) 
D^ 2 = d^ 2 - i A 1 fa - i • (2.25b) 
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Here the anti-Hermitian bundle connection A^ (= — A^) takes its values in the 
four- dimensional Lie algebra u(2) of the two-particle structure group U(2) and 
may therefore be decomposed with respect to some appropriate basis {ti, r 2 ; x-> X.} 
of u(2) as follows 

2 

A, = A%Ta + B, X ~ S; X ■ (2.26) 

0=1 

Clearly a similar decomposition does apply also to the curvature (12. 23ft of 
Afj_, i.e. 

2 

^> = E F V r « + - G>X , (2-27) 

a=l 

where the curvature components (field strengths) are given in terms of the connec- 
tion components (gauge potentials) through 

(2.28a) 
(2.28b) 

A\] B„ (2.28c) 
A\] B; . (2.28d) 

The physical meaning of such a decomposition refers to the subdivision of the 
totality of gauge interactions into the electromagnetic forces (mediated by the 
electromagnetic potentials A a ^ ) and into the exchange forces (mediated by the 
exchange potential B^). Observe here that both wave functions ipi and ip2 be- 
come directly coupled to each other by the exchange potential B^ via the minimal 
coupling (l2^oB - fl2~25bl . 

The interesting point with the exchange forces refers now to the fact, that 
they can act exclusively among identical particles (sharing identical masses and 
charges). Or conversely, for non-identical particles (such as the constituents of our 
positronium system) the exchange forces must vanish (i.e. = 0) so that the 
non-identical particles can feel only the usual electromagnetic forces. The reason 
for this is due to the former requirement of covariant constancy of the total velocity 
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f\ u = v,a\ - v u a\ + i [b,b; - b u b;] 
f\ u = v,a\ - v v a\ - i [b,b; - b v b;] 

Gfiu = V»B V - V V B^ + 1 [A\ - A 2 J B v - i [A\ - 

g% = w,b; - v v b; - i [a\ - a\] b; + i [a\ - 



operator (12.191) and the mass operator M. ( I2.20p . Indeed since these objects are 
of absolute (i.e. non-dynamical) character, they can be thought to be constant over 
space-time and consequently those gauge covariant constancy conditions imply the 
commutativity of those objects with the structure algebra being spanned by the 
above mentioned generators {r a } = {r a ; x, X.}'- 

[M,r a ]=0 (2.29a) 
PT M , r a ] = (2.29b) 
(a = l,...,4) . 

This property of commutativity is safely ensured when both particles have identical 
masses and charges so that the absolute objects become proportional to unity 1( 8 ) 
of the eight-dimensional representation of the Clifford algebra C(l,3) 



t^t u + rvr M = r^r, + r.r^ = 2 9flu • i (8) , (2.30) 

i.e. when the mass and velocity operators adopt the following shape: 

M => M ■ 1 (8) (2.31a) 

KV => = ltl ■ 1 {8) . (2.31b) 



Thus for identical particles, the commutation conditions (I2.29al) - fl2.29bl) are triv- 
ially satisfied and the gauge potential (12.261) can sweep out the whole structure 
algebra u(2), with non-trivial exchange fields B^. 

However when the particles carry opposite charges, the velocity operator IT M 
(12. ip is to be used in place of (= 7 M ©7^); and furthermore, when both particles 
have different masses M a (a = 1,2), the mass operator is no longer proportional 
to unity but adopts the more general shape 

2 

M = iJ2 M aT a - (2.32) 

a=l 
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Naturally for this general situation, the absolute objects ./Vf,F M do not commute 
with all four generators r a of the original structure algebra u(2) which therefore 
must be reduced to a certain subalgebra u'(2) G u(2) so that the commutativity 
requirement does then hold with respect to that subalgebra u'. For the present 
positronium case, this reduction process cuts the original structure group U(2) 
down to its maximal Abelian subgroup U'(2) — U(l) x U(l) which is due to the 
two electromagnetic generators r a (a = 1,2). The latter are obeying the following 
commutation relations: 

[n, r 2 ] = [x, x] = -i (n - r 2 ) 

[n,x] = iX \ T 2, x] = ~iX ( 2 - 33 ) 

[n,x] = -ix [T2,x]=ix- 

According to this reduction process, the bundle connection A^ (12. 26f) and its cur- 
vature (12.271) become projected to the abelian subalgebra u(l) ©u(l), i.e. the 
exchange potential is to be omitted: 

2 

a^J2 A \ r - = A \ T ^ + A V 2 ( 2 - 34a ) 

o=l 
2 

j> =► E F V r « = F V r i + ( 2 - 34b ) 

a=l 

F%, =► V,A\ - V U A\ . (2.34c) 

This situation of missing exchange interaction does apply to the subsequent treat- 
ment of the positronium system consisting of a positron and an electron which, 
it is true, have the same mass M (I2.31al) but are oppositely charged so that the 
velocity operator F M (12.11) is to be applied in place of its identical-particle form 

A similar argument for the reduction of the structure group comes from the 
inspection of the RST currents j afl (a = 1 ... 4) which are quite generally defined 
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in terms of the gauge velocity operators v ajjl through 

j a » = yv*^ (2.35) 

The point with these currents is here that (by explicit differentiation with reference 
;o the RST dynamics) their source equations are found to be of the following form 

\ 

D»j a „ = -^-J?[Mc 2 ,T a ]q . (2.36) 

But since these RST currents act as the sources of the gauge fields A a (see the 
non-Abelian Maxwell equations below), the currents j afl must themselves have 
vanishing source from consistency reasons: 

D»j ap = . (2.37) 

However, this requirement can be satisfied again by reducing the structure group in 
such a way that the commutator f)2.29ap vanishes and the desired gauge continuity 
equation ( 12.371) becomes then implied by the general relation (12.361) . 

Thus, the suppression of exchange forces among non-identical particles is seen 
to be supported by the intrinsic logic of RST, whereas in the conventional theory 
the emergence/suppression of exchange effects for (non-)identical particles is a 
purely kinematical effect being induced by postulating the (anti)symmetrization 
of wave functions or (anti) commutation of field operators, resp. 

E. Gauge Field Equations 

The RST dynamics presented up to now is not yet a closed system because 
one finally has to specify a field equation for the gauge potential A^. The choice 
of such a field equation must necessarily undergo certain constraints because it 
must be compatible with the already specified RST dynamics, especially with 
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the conservation laws. Our choice is the non-Abelian Maxwell equation for the 
curvature (I2.23P 

OT^ = -47Tia s X. (2.38) 

Here as is the fine structure constant (e 2 /hc); and the current operator may 
be decomposed with respect to the chosen Lie algebra basis {r a } as 

c7 A4 = *j>« (2.39) 

so that the abstract equation ( 12. 38ft reads in component form 

D^F^ v = Ana s f v . (2.40) 

As a brief consistency check of the chosen gauge field dynamics, convince yourself of 
the fact that those continuity equations ( 12. 37ft are an implication of the Maxwellian 
equations ( 12. 38ft . or (12.401) . resp. To this end, simply apply the bundle identity for 
the curvature 

V^T^ = (2.41) 
to the Maxwell equations ( 12. 38ft and thus find the desired source equations 

V 1 3^ = D^j a ^ = . (2.42) 

It is true, these source equations (12.421) for the Maxwell currents j a are not 
strictly identical to the corresponding source equations ( 12. 3 7ft for the RST currents 
jan (12.3511 . But it is possible to introduce a fibre metric K a p in the Lie algebra bun- 
e, which is invariant under the structure group and is also covariantly constant 
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2fl 

D„K aP = . (2.43) 

This, namely, provides one with the possibility to identify the Maxwell and RST 
currents as contra- and covariant versions of each other: 

j% = K^jp, (2.44a) 
] a » = K aP j% . (2.44b) 
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Clearly through this identification, both forms of continuity equations (12.371) and 
(I2.42p become now equivalent and this does support the desired compatibility of 
the matter and gauge field dynamics. 

For the sake of generality, one may first wish to include the exchange interactions 
(i.e. putting ^ 0) so that the Maxwell equations (I2.40p read explicitly 

WF\ U + i [B»G%, - B^G, U ] = Arra s] \ (2.45a) 

VF%, - i [B»G% - ffVG,.,] = 4vra s j\ (2.45b) 

VG^ + i [A\ - A\] G\ - i [F\ v - F\ v ] B» = 4vra s g v (2.45c) 

^G; u - i [A\ - A\] G*» + i [F\ u - F\ v ] B^ = Arra s £ , (2.45d) 

where the exchange currents j 3 ^, j 4 ^, (I2.39H are denoted by 

3\ = 9, (2.46a) 
j\ = -g; • (2.46b) 

However, since the positronium system consists of two non-identical particles which 
are not able to feel the exchange force, one has to drop the exchange potential 
and thus the gauge field system (I2.45al) - (12.45dl) simplifies to 

V»F\ v = A<Kasj\ (2.47a) 
V"F%,= Arras ]\ • (2.47b) 

From this reduction process it is clearly seen that the omission of the exchange 
force implies the simplification to an Abelian and linear theory referring to the 
electromagnetic interactions alone. 

Finally, in order to close the whole dynamical system, one has to specify the 
link of the Maxwell currents (a = 1,2) to the one-particle wave functions ip a . 
This may be done by first defining the Dirac currents k aii in the usual way through 
(a = 1,2) 

= i'al^a , (2.48) 
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and then the RST currents j afl ( I2.35P are found to be related to these Dirac currents 
through 



h» = km (2.49a) 
hv. = ■ (2.49b) 

Thus, the last step for constructing the Maxwell currents j a ^ (I2.44aj) must 
consist in specifying the fibre sufanetric K ab or K ab , resp. In general, this object 
depends upon one parameter u (i.e. the self-interaction parameter): K ab = K ab (u). 
However for the present purposes, we want to neglect the self-energy effects which 
implies u = 0; and for this special value of u the two-dimensional fibre submetric 



becomes |z] 



= {K ab (0)} = \ \ ^ \ (2-50) 



so that the desired two Maxwell currents j a ^ read in terms of the wave functions 

3 X y, = -hi* = k m = 4>il^i (2.51a) 
ff, = -3i,m = ~k 2f , = -$21^2 ■ (2.51b) 

Thus the first Maxwell current j 1 ^ correctly refers to a positively charged particle 
(positron) and the second Maxwell current j 2 M refers to the negatively charged 
electron. 

In this way, the RST dynamics for a two-particle system becomes closed in 
a consistent way: the wave function of the first /second particle generates the 
first/second Dirac current fei^/A^ according to (12.481) ; then these Dirac currents 
^a/i gi ve r i se to the Maxwell currents according to fl2.51al) - fl2.51bl) . Next, the 
Maxwell currents act as the sources for the generation of the gauge potentials 
A a ^ as the solutions of the (Abelian) Maxwell equations (I2.47ap -( l2.47bj) : and finally 



the first/second gauge potential A 1 /A 2 enters the Dirac equation (12.181) for the 
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second/first particle via the gauge-covariant derivatives fl2.25ap - fl2.25bj) . Clearly 



such a simple cross relation is to be expected from logical reasons, but it must be 
somewhat modified when the self-interactions and exchange effects are included 

ik 

III. MASS EIGENVALUE EQUATIONS 

In a stationary bound state, both positronium constituents (i.e. positron and 
electron) will be described by wave functions ip a (r,t) satisfying the usual space- 
time factorization 

r M*c 2 i 

ipi(f, t) = exp i — ^ — t ■ ipAr) (3.1a) 
L n J 

r M*c 2 i 

fair, t) = exp -i — — t -fatr) • (3.1b) 
L n J 

Observe here that the time factors of both constituents are inverse to each other 
which is due to the fact that the first particle (a = 1, positron) obeys the Dirac 
equation (12.181) with a negative mass term 

ih^D^t = -Map! , (3.2) 

whereas the electron obeys the conventional Dirac equation 

ih-fDnfo = Mc*{j 2 . (3.3) 

Clearly, this difference is due to the former choice of the total velocity operator 
F M (12.1 p for oppositely charged particles. 

A. Positronium Configurations 

The mass eigenvalue M* in (I3.1al) - (13.1bp is common to both particles and must 
be determined from the coupled Dirac equations ( I3.2p -( l3~3j) by use of the stationary 
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ansatz (I3.lap - fl3.lbp . Here it is very advantageous to split up the Dirac four-spinors 
?/> a (r) into the direct sum of two- component Pauli spinors ^^>±(r), i.e. one puts 

^1 = (a WW© (a V_(f) , (3.4) 

and then the mass eigenvalue equation for the first particle (positron) reads as 



follows 



i (a> V)<V(fO + (2) ^o ■ (1 Wr) - M • *< \ ± (f) = ±M + M * c . (D^ (f) , (3.5) 

and similarly for the second particle (electron) 

i (a • V)( 2 V±(r) + {1 U • < 2 WrO - ^1 • ^±(0 = -^^c • ( 2 V T (f) . (3.6) 

ft 

Now the crucial point with the positronium system is that its energy spectrum 
is experimentally found to be essentially a one-particle spectrum [181 ]. the energy 
levels of which are to be subdivided in RST into two subsets: ortho-positronium 
and para-positronium. These two classes may be characterized by the relative 
orientation of the magnetic fields H a (r) (a = 1, 2) being emitted by both particles: 



ortho-positronium : Hi(r) = H2{r) = Hb(r), Ai(r) = A 2 (r) = A&(r) (3.7a) 
para-positronium : Hi(f) = —H2(r) = H p (f), Ai(r) = — Azif) = A p (f) . 

(3.7b) 

In contrast to this, the conventional classification refers to the spin orientations 
and thus works rather in terms of the singlet ( 1 So) and triplet ( 3 Si) states 

Since, however, the magnetic fields are after all generated by the Dirac three- 
currents k a (r) = {^k j (r)}, cf. (gHHD 

k a {r) =$ a {r)l^a{r) , (3-8) 
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and since both particles carry opposite charges one will expect the following ar- 
rangement of the three-currents: 

ortho-positronium : ki(f) = —k2(f) = k p {r) (3.9a) 
para-positronium : k\{f) = k2(r) = kb(r) . (3.9b) 

Actually in view of such an arrangement, the corresponding Maxwell three- currents 
J a {r) fl23Ta) - fl23Tbj) become 



ortho-positronium : ji(r) = j 2 (r) = = k p (r) (3.10a) 
para-positronium : = —J2{t) =f j P {r) = kj,{r) . (3.10b) 

Such an arrangement of the Maxwell currents is immediately plausible, because 
this kind of currents is responsible for the generation of the magnetic fields H a (r) 
according to the (Abelian) Maxwell equations (I2.47aj) -( l2.47bj) . the space part of 
which reads in three-vector notation 

V x H a = Anas Ja (r ) • (3.11) 

Thus the (anti)parallelity of the Maxwell currents j a (r) (I3.10al) -( l3.10bl) actually 
implies the (anti)parallelity of the mag netic fields H a (r) fl3T7aD - fl3~7bj) 



On the other hand, the situation with the electric fields E a (r) is less com- 
plicated: since both particles are oppositely charged, the Dirac densities ^k (r) 
(a = 1, 2) must be identical 

ortho & para-positronium : ^Afo(r) = ^k (f) = ^ b \ (r) , (3-12) 

because in this case the Maxwell densities ^ a 'jo(r) f l2.51al) -( [2.51b[) as the true charge 
densities will differ in sign as required for oppositely charged particles: 

ortho & para-positronium : jo(r) = jo(^) =f jo(r) = ko(r) . (3.13) 

As a consequence, the electric fields E a (r) will also differ in sign 

ortho & para-positronium : Ei(r) = —E 2 (r) = E p (r) . (3-14) 
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Furthermore, the missing of exchange interactions lets the electric fields E a (f) 
appear as pure gradient fields, cf. (I2.34cl) 

4(r) = -VH4 (r) , (3.15) 

with the electrostatic potentials ^Ao(r) to be identified in the following way 

ortho Sz para-positronium : ^Ao(r) = —^Ao(r) = ^Ao(r) . (3.16) 

Indeed, this electrostatic arrangement matches again the (time-component of) 
Maxwell's equations fl2^?aT) - (l2.47bh 

V-4 = -47ra s (a) j'o(r) > ( 3 - 17 ) 
i.e. under use of the electrostatic identifications ( l3.13l) -( l3TT6l) : 

A (p U (f) = -47ra s (p) j (r) . (3.18) 
B. One-Particle Equations 

This specific way of identifications of densities for both positronium config- 
urations must now necessarily imply a corresponding identification of the Pauli 
spinors ( a «p±(r) (13 Ah since the charge and current densities are ultimately built 
up by the wave functions if> a (f). The result of such a spinor identification is the 
contraction of the four mass eigenvalue equations (I3.5l) - (l3.6p to only two equations 
which then formally describe only one particle. This is the RST analogue to the 
conventional energy eigenvalue equation being written in terms of the relative co- 
ordinate f (= r*i — F2) after the center-of-mass motion R (— fi + r*2) has been 



separated-off [10]. 



Turning first to the simpler case of the Dirac densities ^^(r), one finds these 
objects emerging in terms of the Pauli spinors *-°V± (13.4j) as 

{a) k (r) = ipaloA = (a H(r) • W<p + (r) + (<V_(r) • ^{r) . (3.19) 
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The somewhat more complicated case of the Dirac three-currents k a (f) (j3.8p looks 
as follows 

K(r) = Mr)lMr) = {a V + (r)a^_(f) + (<V_(f) a^ + {r) , (3.20) 

where 7 denotes the space part of the Dirac matrices 7^ (standard form) and 
<? — {o~ x , &y, o~ z ) refers to the well-known Pauli matrices. 

Now, in order to let the set of four two-particle equations (I3.5p - fl3.6l) collapse 
to the desired set of two one-particle equations under observation of the required 
density and current identifications, one expresses the second Pauli spinors ^<f±(r) 
in terms of the first Pauli spinors ^\p±{r) (4= ^±(^0) i n the following way 

ortho-positronium : ( - 2 ^ + (r) = Q < - 1 ^ + (r) = Q y?+(r) (3.21a) 

%_(r) = = -Q<p-(r) (3.21b) 

para-positronium : ^ 2 V±(r) = Cl^ip±{f) = Qip±(r) . (3.21c) 

Here, the operator Q acts in the two-dimensional Pauli-spinor space and is to be 
determined from the identifications of densities and currents. Its first property 
is found from the identification f)3.12l) of the Dirac densities ^ a ^ (r) f)3.19p . which 
holds for the ortho- and para-configurations. Evidently, this density identification 
reads explicitly 

^(r)(fi f fi)^+(f) + y? t _(f)(7i t ^_(f) = <pl(r)tp + (r) + ipUr)V-(r) (3.22) 

and thus is surely obeyed by adopting the operator Cl to be unitary ($V = Q~ r ) 

&(l = l. (3.23) 

The next property of the operator Cl may be deduced from the current iden- 
tifications (I3.9ap - fl3.9bp which by use of equation ( 13.201) explicitly read for both 
configurations 

= ^ + {r) a(p-(r) + (pL(r) a<p + (r) . (3.24) 
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But this requirement can easily be satisfied by adopting the following additional 
property of fl: 

&(k-a)fl = k-a (3.25) 
(ortho & para-positronium) 

where the unit three- vector k (\\k\\ = 1) points in the direction of the corresponding 
three-current, i.e. one has (cf. (^a] - (j3T9bl ) 

ortho-positronium : k p (r) = k p (r) k(f) (3.26a) 
para-positronium : kbif) = kbif) k(f) . (3.26b) 

But clearly, it is no problem to find unitary solutions (13.231) for the equation (13. 25j) . 
namely the following two 

ortho-positronium : Cl — i ■ 1 , (3.27a) 
para-positronium : Cl = i k* a . (3.27b) 

We shall readily elaborate the reason for this peculiar arrangement of the ortho- 
and para-cases, resp. 

But now that the operator fl is known, one can turn to the main problem, 
namely the reduction of the two-particle equations (I3.5l) -( l3~6i) to one-particle equa- 
tions by means of the preceding identifications. Considering first the ortho-case, 
one finds the mass eigenvalue equations (13.61) of the second particle (electron) 
adopting the following form through application of the relations fl3.21ap - fl3.21bj) : 

i(a*V)£lip+(r) -^A -Cl<f4r) -A b *a£lp + (r) = M * ^ M c- (l^_{r) 

(3.28a) 

-i(a-V)n V _{r) + ^Av-n V+ {r) + A h -an V _{r) = - *~ c-ttip + (r) , 
v ' n 

(3.28b) 
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whereas the ortho-identifications recast the eigenvalue equations (13.51) of the first 
particle (positron) to the following form: 

z(a-VWf) - {p) A -<fJr) - A b -c<p + (r) = M * ^ M c-^_(r) (3.29a) 

' 1 h 

i(a*V)v-{r) - {p) A ■ <p+(r) - A b ' a ipJr) = M * ~ M - c • <p + (r) (3.29b) 

( ortho-positronium ) . 

Hence the identification of both sets of equations ( I3.28al) - (13.28bl) and (I3.29al) - 
( 13.29bl) yields a further condition upon the operator Cl, namely 

fi f (a- V)Cl<p±(r) = a- V (pjf) (3.30a) 
&(A b -a) Cl = (A h 'B) (3.30b) 
( ortho-positronium ) . 

Evidently, these two requirements for ortho-positronium do admit the solution 
(I3.27aj) for the operator Cl with no further constraints which is in contrast to the 
corresponding situation with the para-configuration. 

The treatment of the para-configuration runs in a quite analogous way: first, 
apply the para- transformation (13. 2 left again to the eigenvalue equations (13.61) for 
the electron and find in this case 

z(a-V)fi^ + (f) + (p U -fi^(f) - (A p *a)(l<p + {r) = - c-£lip-{r) 

(3.31a) 

i(a-V)n<p-(r) + {p) A -n(p + (r)- (A p - a)O^-(r) = - — ~ M c ■ ^y? + (r) , 

(3.31b) 

whereas the eigenvalue equations H 3 . 5 f) of the positron appear under the action of 
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the para-identifications as 

V)^ + (r) - ( p l4 • <fi-(r) + (V5> + (r) = M * + M C -y?_(r) (3.32a) 

z(c? • V)v5-(r) - (p U • p+(f) + (4 • a)^_(r) = M * ~ M C • p+(f) (3.32b) 

( para-positronium ) . 



Therefore, if both sets of eigenvalue equations fl3.31ap - fl3.31bp and f)3.32ap - fl3.32bp 



are required again to be identical, one is led in this case to the following condition 



&(a- V)n<p±(r) = -(<?• V)p±(r) (3.33a) 
A t (^-?)fi= (Ap'ff) (3.33b) 
( para-positronium ) . 

Observe here that the first condition f)3.33al) does not admit the first solution 
f)3.27al) for Cl, which therefore can apply only to the ortho-configurations. But 
the requirement f)3.33al) does admit the second solution (13.27bl) ; and in this case 
the relation (I3.33al) presents some restrictive condition upon the form of the Pauli 
spinors <fi±{f) (i.e. double valuedness, see below). Furthermore, the second con- 
dition (13.33bp with the operator Q being given by (13.27bp requires the vector 
potential A p for the para-configurations to be (anti)parallel to the Dirac current 

A p = A p k . (3.34) 

Thus, summarizing the situation with the one-particle eigenvalue equations, 
one arrives at two forms thereof: namely, the ortho-form f l3.29ap -f l3.29bp and 
the para-form fl3.32ap -f l3.32bp . Furthermore, the ortho-form is found to admit 
a richer spectrum of mass eigensolutions; the reason is that the shape of the ortho- 
eigensolutions (p±(r) is not constrained by the ortho-condition f)3.30ap . Naturally, 
further insight into this peculiarity of the positronium eigenvalue problem will now 
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be gained by a closer inspection of just those restrictive conditions (I3.33ap - fl3.33bj) 
upon the form of the para-eigensolutions f±(r). 



C. Para-Positronium 

The task is now to elaborate the specific form of the Pauli spinors which is 
admitted by the condition f l3.33al) with the operator Q being given by equation 
(13.27bp . For the sake of simplicity, we assume a cylindrical symmetry so that both 
Dirac currents k b / p f)3.9ap and (13.9bp encircle the z-axis, where the unit vector k 



coincides with the azimuthal basis vector of an orthonormal triad {e r , e&, e^} 
being due to the spherical polar coordinates {r, $,0}, i.e. we put for the vector 
potentials 

A b/p (r) = ( - b ^A (j> (r ) <d)e (t>) (3.35) 
and similarly for the Dirac currents 

k b/p (r) = ^%(r,$)e ( p, (3.36) 

cf. fl3.26ap - fl3.26bp . By this assumption of cylindrical symmetry (for a toroidal 



symmetry see ref. [l5|), the para- condition fl3.33al) adopts the following anticom- 
mutator form 

{(a>V),e <P 'a}<p ± (r)=0 (3.37) 
( para-positronium ) . 

But here one can show by means of some elementary mathematics, that this 
condition upon the Pauli spinors (f±(f) can be recast to the following form 

{i a • (V x + V • e +2(e • V)}p±(r) = (3.38) 
( para-positronium ) , 

where of course the divergence of the azimuthal basis vector vanishes 

V-e^O, (3.39) 
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and the curl of es is given by 



V x e = . (3.40) 

rsmw 



Thus the restriction (I3.38P upon the Pauli spinors <p±(r) adopts its final form as 

lMr) = (3.41) 
( para-positronium ) , 

where the z-component of the total angular momentum operator J(= S + L) is 
defined as usual through 

Consequently, the physical content of the para-restriction (I3.33al) is nothing 
else than the fact that the para-eigensolutions must have vanishing ^-component 
of total angular momentum. Observe that this requirement refers to a one-particle 
spinor state ip(r) = <p+(r) © <fi-{f) which is the direct sum of both Pauli spinors 
(p+(r) and (p_(r). Such exotic states do not occur in the conventional theory of a 
single spin- 1 Dirac particle just on behalf of the required uniqueness of the wave 
function; but in the present positronium description, these exotic states must nec- 
essarily emerge as a consequence of the physical equivalence of both constituents, 
cf. the preceeding identification process leading to (13.33aj) . 

It is also interesting to remark that para-positronium is somewhat exceptional 
also in other respects: since the total charge density j is zero according to equation 

{EH 

Jo = (1) Jo + (2) Jo = , (3.43) 
as well as the total electric potential Aq, cf. (13.161) 

A = WA + ^A = , (3.44) 

the total electric field strength E must also be zero: 

E(r) = E^r) + E 2 (f) = -VA (r) = . (3.45) 
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A(r) 



Similar conclusions do hold also for the total magnetic objects, i.e. 

= Ji(r)+J 2 (r) =0 (3.46a) 
= li(r) +A 2 (r) =0 (3.46b) 
#(f) = #i(r) + H 2 (r) = , (3.46c) 

see equation ( 13.7bl) and ( 13.10bl) . On the other hand, since the composite RST 
systems do interact with external sources just via their total objects (see the 



15|), the para-positronium 



discussion of this point in the preceding papers 
appears to be reluctant to undergo interactions with the outside world. 

Concerning now the construction of appropriate Pauli spinors if(f) obeying 

the angular-momentum condition (I3.4ip . one first constructs four basis spinors 

r..(+) ,,(-).,.(+) ,,H\. 

\U) ,CU ,LU 1 ,LU 1 }. 



UJ, 



(+) 



(-) 



-i(f>/2 z-2'2 
S 



A6/2 ri~2 



Co 



U! 



(+) 



(-) 



l l 

e -l0/2 £2 < 2 



J4>/2 /-2' 2 



C 2 i 



(3.47a) 
(3.47b) 
(3.47c) 
(3.47d) 



T < ; (±) 



which are required to have vanishing z-component of total angular momentum J z : 

J z uj^ ] = 0. (3.48) 
This goal can be achieved by taking for the £-basis the conventional choice [3] 

•/Vr lr.li.1 ■ 1K7' (3.49a) 

./ c j, r n,h er (3.49b) 

t C j f = 1(1 + 1) CT (3.49c) 
I 2 = h 2 s(s + 1) (T , (3-49d) 

with the following spin-| combination of quantum numbers for the positronium 

groundstate 

s = j = \; m = ±\; 1 = 0,1. (3.50) 
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And then one decomposes the Pauli spinors <p±(r) with respect to this u;-basis as 

cp + (f) = R + ( r , ■&) ■ 4 +) + S+(r, ■&) ■ 4 _) (3.51a) 
ip_{f) = -iR_(r,&) -uj{ +) -iS-(r,ti) ■ . (3.51b) 

Finally, one inserts this form of the Pauli spinors into the eigenvalue equations 
fl3.32al) - (13.32bl) for para-positronium in order to deduce the corresponding eigen- 
value equations for the wave amplitudes R±(r, i?) and S±(r, i?) as [9(: 

3R+ ldS + (vU ~ (v) . r . a ~ ~ . M + M* ~ 

— - + -- ^± - (P U ■ A_ + sin ■ i?+ - cos t? ■ S + ] = c ■ R. 

or r ov n 

(3.52a) 

ijjrkA _ i_as_ + k _ ms§ = m-m 1c 

r or r ov n 

(3.52b) 

^ - - (p U • S_ - H^sin^ • S + + costf • = ^i±Ml c . s_ 

or r ov n 

(3.52c) 

UMA + 1_BR_ + w . s + w # . s _ + cos , . ^ = M_M lc . ^ 
r or r ov n 

(3.52d) 



Here, the original wave amplitudes -R±(r, i?) and S±(r, $) f)3.51ap - fl3.51bp have been 
slightly changed to 

R ± (r, tf) = Vr sin ■ R±(r, 0) (3.53a) 



S±(r, i?) = Vr sintf • 5±(r, ti) ; (3.53b) 

and furthermore, since the vector potential A p (13.7bl) must always be (anti)parallel 
to the unit vector k = (13.341) . we have put also for its azimuthal component: 

A p = V>A& . (3.54) 

Notice also that the use of the cu-basis f l3.47al) - (13.47dl) induces a certain non- 
uniqueness of the Pauli spinors <f±(r) fl3.51ap - fl3.51bp ; indeed it is obvious that 
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they are unique over space-time only up to sign: 

¥>±M,0 + 27r) = -p ± M,0) . (3.55) 

However, it is also obvious that the physical densities, e.g. charge density ^k (r) 
(I3.19P and current k a (r) f)3.20p . are nevertheless unique! 

D. Ortho-Positronium 

Naturally, the ortho-configuration must appear somewhat simpler because the 
corresponding operator Cl ( 13.27aj) does trivially obey the ortho-constraint (I3.30aj) 
and therefore does not impose any restrictive condition upon the Pauli spinors 
l P±{f)i as is the case with the para-configuration (I3.33aj) . Nevertheless, it is rea- 
sonable to assume that also the stationary ortho-configurations are associated with 
a well-defined z-component J z of total angular-momentum; albeit this component 
is not necessarily zero as for the para-configurations. But for the special case 
of the groundstate it is self-suggestive to adopt zero angular-momentum also for 
ortho-positronium. This means that we can resort to the above mentioned w-basis 
(E3IaD-(j33zaj) and thus the para-ansatz ( 13.51ap -l l3.51b)) may be applied also for 
the ortho-configurations! Consequently, one inserts this ansatz into the ortho- 
system ( I3.29al) - (13.29bj) and then arrives at the following eigenvalue system for the 
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two wave amplitudes R±(r,i9) and S±(r, 

dR+ ldS+ (vU ~ m . . . a ~ ~ . M + M* ~ 

+ ~- {P) A ■ i?_ - sin v ■ R+ - costf • S+] = c ■ R_ 

or r ov n 

(3.56a) 

1 0(rJL) _ 1 d§_ + w ^ + + w ^ _ co ^_ - = M-M, c ^ + 

r or r ov ' a 

(3.56b) 

-t^ ;ttt - • 5- + % sin i? • S + + cos i? • R+] = c • S_ 

Or r ov n 

(3.56c) 

1 d(rSS) 1 OR. Ma ~ m . r . ~ - . M - M* ~ 

— + — sir + A o • 5+ - ( % sin i? ■ S_ + cos i? ■ R_] = c • 5+ . 

r Or r Ov n 

(3.56d) 



Evidently, this ortho-system differs from its para-counterpart (13.52aj) - (13.52d|) 



merely in the sign of the magnetic term (^A^ — ► —^A^); and if the magnetic 
(i.e. spin-spin) interactions are neglected, both systems collapse into the same 
form ("electrostatic approximation"). This is reasonable because the ortho-para 
dichotomy of positronium is just due to the different magnetic interactions and 
therefore must disappear when the latter type of interaction is neglected. This 
circumstance will readily be exploited for a lowest-order approximation. 



E. Approximations 

It should be obvious that both eigenvalue problems (j3.52aj) - (13.52d[) and fl3.56al) - 
( I3.56dl) are too complicated in order to find exact solutions (which is mostly the 
case with modern gauge field theories). Therefore it is important to look for simple 
approximative solutions which, however, are still well-suited to demonstrate the 
essential physical effects. Surely, the binding phenomenon can be understood al- 
ready by reference to the electric interactions alone, which mostly in atomic physics 
are much stronger than their magnetic counterparts. Therefore one may expect 
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to obtain a relatively good approximation to the binding energy by neglecting the 
magnetic terms in the eigenvalue equations, with preservation of their relativistic 
form. Furthermore it seems reasonable to assume that magnetism is closely re- 
lated to a non-trivial angular momentum, which however is known to be in conflict 
with the spherical symmetry of the field configuration. Therefore, if one wishes to 
consider first a "spherically symmetric approximation" of the positronium ground- 
state, one will skip the ^-dependence of the wave amplitudes R±, S± together with 
the magnetic (spin-spin) interactions. By these approximative assumptions both 
eigenvalue systems (I3.52al) - (13.52dl) and (13.56aj) - (13.56dj) coincide to the same form: 



alR + {r) 

dr 
ri?_(r)) 
dr 

dS + (r) 
dr 





r) • R„\ 


v r)= h c 


R. 


(r) 


(3.57a) 


+ {P) A ( 


r)-R + i 


K r)= h c 




(r) 


(3.57b) 


- (p U i 


(r) ■ S_i 


, , M + M* 
v r)= h c 




(r) 


(3.57c) 


+ 


(r) ■ S + > 


,r)= h c 


■S + 


(r). 


(3.57d) 



r dr 

The striking feature of this simplification evidently refers to the circumstance 
that the spin-up fields R±(r) become decoupled from the spin-down fields S±(r); 
and therefore it is sufficient to consider the spin-up configurations R±(r) alone, 
which are solutions of the coupled pair of equations (I3.57al) -( [3.57b[) . Clearly, the 
spin-down solutions S±(r) of (13.57c| - (13.57dj) are equally well possible and are iden- 
tical (up to spin-orientation) with the spin-up solutions. Consequently we can 
concentrate exclusively upon the spin- up solutions R±(r). 

The next problem refers to the electrostatic potential ( p l4 (r) which has also 
been assumed to be spherically symmetric in the simplified system (I3.57a!) -( l3.57dl) . 
Properly speaking, this potential must be determined from its own Poisson equa- 
tion, cf. (13.181) . because it is generated by the charge distribution ( |3,13[) of 
the other particle. However, it is a well-known fact in classical electrodynamics 
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that the potential ^Ao(r) as a solution of the Poisson equation (13. 18ft must be of 
the following form 

^„(r)=a s /dV'^)^^, (3.58) 
J \\r — r'\ \ (r-+oo) r 

which in the asymptotic region (r — >oc) adopts the Coulomb form (~ r" 1 ). This 
fact may now be exploited in order to briefly demonstrate that the residual eigen- 
value equations 

<!Ml-?±. k _ (r) = E^k c . k _ {r) (3.59a) 
dr r n 

l d(rR.(r)) g s M - M. ~ 

— ; - H ■ R+(r) = c-R+lr) (3.59b) 

r dr r n 

do not agree with the well-known solutions of the conventional Coulomb force 
problem, despite the use of the asymptotic Coulomb potential (I3.58p . Indeed, 
trying here the usual ansatz for the wave amplitudes R±(r) 

R ± (r) = N ± ■ r u exp[~— } (3.60) 

yields by insertion into the simplified eigenvalue equations (I3.59aj) - ( 13.59bl) for the 
scale parameter r*, power u, mass eigenvalue M* and normalization constants N± 
the following results: 

1 h 2 

r* = -a B = 5 • • • Bohr radius (3.61a) 

2 ZMe z 

v = -\{l-^{l-{2a s f) (3.61b) 



M* = MV1 - (2a s ) 2 (3.61c) 



N- 1 - a/1 - (2a s 2 

^_ = Y 3.61d 

But these results do not agree with the corresponding ones of the conven- 

n 

tional Coulomb force problem (see, e.g., ref. [20J), which may be seen even more 
clearly through passing over to the non-relativistic approximation. For this ap- 
proximation, the scale parameter r* ( I3.61al) remains unchanged whereas the power 
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v (I3.61bj) is put to zero on account of the smallness of the fine structure constant 
«s = e 2 /he {pa 1/137). Furthermore, the negative Pauli component i?_ is put to 
zero for the same reason (see the ratio ( I3.61d[) of the normalization constants N±), 
and finally the mass eigenvalue M* is expanded with respect to as 

Me 4 

M* = M(l - 2a 2 s ) = M — 2 — - . (3.62) 

nrcr 

Defining here the non-relativistic eigenvalue E$ through 

E s = M,c 2 - Me 2 (3.63) 

yields finally 

Me 4 e 2 
E s = -2— = -2—. 3.64 

nr a-Q 

Thus the groundstate binding energy E B (= —Eg) of a point particle of mass M 
and charge e in the Coulomb potential (13.581) is found here to be four times the 
conventional hydrogen result 

o2 , „2 



E f 



= 2 — = 4 • Eb 

rst a B 



(For the conventional treatment of the hydrogen atom see, e.g., ref. 20|). The 
non-conventional RST result (I3.65P does not mean that RST is wrong, but simply 
says that RST treats the internal motion of a compound system in a way different 
from to the external motion of a pointlike particle! 

In order to complete this RST picture of the internal motion of the positronium 
constituents, one wishes to see the non-relativistic approximation (R{r), say) of 
the wave function R±{r) (13.601) 



R+{r) =^ iVexp 



= R(r) (3.66a) 



r " 
r*. 

R_{r)^0. (3.66b) 



Naturally, one expects that this non-relativistic form of the solution will satisfy 
some non-relativistic wave equation (for the internal Coulomb force problem); but 
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this wave equation cannot coincide with the conventional Schrodinger equation 
because of the non-conventional eigenvalue Eg (I3.64p . Rather, the desired non- 
relativistic form of the original RST eigenvalue equations (I3.59al) -f f3.59bl) . with 



arbitrary electric potential Ao(r), is obtained by approximating the first one (I3.59ap 
through 

Ur) - 4r ■ ^ (3-67) 
v ; 2Mc dr v ; 



and substituting this into the second one f!3.59bl) in order to find 

h 2 ( d 2 id 



2M\dr 2 



(dF? + r^) k{r) ~ hciP)Mr) = Es ' R{r) ■ (3 - 68) 



Clearly, the non-relativistic form fl3.66aj) of the relativistic solution (I3.60p - fl3.61dj) 



together with the Schrodinger eigenvalue Es (13.641) does obey this non-relativistic 
wave equation (13.681) with being identified as the Coulomb potential ( 13. 58ft ; 



but nevertheless this non-relativistic eigenvalue equation (I3.68P is not the ordinary 
Schrodinger equation for the Coulomb force problem! The reason is that it is not 
the Laplacean A(r) = V 2 , in spherical polar coordinates 

dr 2 r dr ' 

which acts upon the wave function R(r). But rather in equation (13.681) . there 
emerges the Laplacean written in cylindrical coordinates, so that the variable r 
changes its geometric meaning from a spherical polar coordinate to a cylindrical 
coordinate. This change then results in the unconventional eigenvalue Es (13.641) . 
in contrast to its conventional counterpart (13.651) . 

Summarizing, the requirement of physical equivalence of the electron and 
positron has lead us to the introduction of those exotic states with vanishing 
z-component J z (13.411) of the total angular momentum; and these exotic states 
were then revealed to obey a wave equation whose non-relativistic approximation 
does not exactly agree with the conventional Schrodinger form. Therefore, when 
the exotic states are used, the non-relativistic limit of RST does not agree with the 
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Hartree-Fock approach which is based upon the Schrodinger type of wave equation. 
This opens the possibility for RST to overcome the deficiencies of the HF approach 
in the sense that the RST predictions may come closer to both the experimental 
data and the conventional Schrodinger predictions than it is possible for the HF 
approach. (For the emergence of the HF approach as the non-relativistic limit of 
RST by using the conventional states, see ref. 13|). 

However there is a second peculiarity of RST by which its non-relativistic limit 
differs from the conventional HF approach: this refers to the specific form of the 
electric potential ^Aq(t), by means of which the positronium constituents are 
attracting each other. For our preceding heuristic demonstrations in connection 
with the approximative wave equations of relativistic ( I3.59al) -( 13.59bl) and non- 
relativistic form (13.681) we preferred to work with the (exact) Coulomb potential 
(( p l4 ^S-; see equation (13.581) ) because this is used also in the conventional 
Schrodinger theory, cf. the non-relativistic positronium Hamiltonian H ( I2.22p . 
However the exact RST potential &'Aq ( 13.58!) surely differs from its asymptotic 
Coulomb form, but not in that way as is the case in the HF approach. Whereas, in 
the latter theory, the singularity of the Coulomb potential (for r — > 0) is completely 
regularized, the RST modification of the Coulomb potential remains singular, but 
in a less pathological manner. Since the specific form of the interaction potential 
( p l4 influences the mass/energy eigenvalues not less than the general form of the 
eigenvalue equations themselves, it is now necessary to inspect thoroughly that 
new form of the electrostatic potential. 



IV. POISSON EQUATIONS 

As the preceding discussion is mainly concentrated upon the eigenvalue equa- 
tions for the matter fields, it must be complemented by a similar investigation 
of the gauge fields, i.e. the electrostatic potential ^l4o(r). This object medi- 
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ates the electric interaction of the positron and electron and thereby dominates 
the magnetic interaction (which is mediated by the azimuthal component ^^^A^ 
of the vector potential Af,/ P (r) (13.351) ). Especially the electric potential is re- 
sponsible for the occurence of binding while the magnetic interactions do at most 
contribute certain small corrections to the binding energy. We therefore neglect 
now the magnetic forces altogether (i.e. putting A p /b to zero) and thus work with 
the truncated eigenvalue systems (13.57aj) - (13.57d|) or A3. 59aj) - A3. 59b|) . resp. These 
equations demonstrate clearly the specific way in which both particles do orga- 
nize their mutual coupling: first, the wave function ip a of either particle builds 
up the corresponding Dirac current k ap , (12.481) . then its Maxwell counterpart j a ^ 
f l23TaT) -f l23Tbl acts as the source of the electromagnetic vector potentials A a 
via the Maxwell equations (I2.47al) - (12.47bl) . and finally these four- vector potentials 
influence the wave function tp a of either particle via the mechanism of minimal 
coupling A2.25aH2.25bh . For the present stationary situation, the Maxwell equa- 
tions (I2.47al) -( [2.47bl) are transcribed to the electrostatic Poisson equation ( 13 . 1 8[) 
for the electric interaction potential ^Ao(r). It is true, this interaction mechanism 
is of the conventional electromagnetic type; but the new aspect is now that the 
potential ^'Aq{t) does adopt a very unusual form when it is generated by those 
exotic one-particle states with vanishing z-component J z of their total angular 
momentum J, see the discussion of equation (13.411) . 



A. Non-relativistic Approximation 

In order to elaborate this point in some detail, one will not be satisfied with the 
asymptotic Coulomb form of the solution ( 13. 58ft to the Poisson equation (13. 181) . 
but one will try to compute that integral (13.581) more rigorously. To this end, 
one makes use of the fact that the Maxwell charge density ^jo(r) (13.131) coincides 
with the Dirac density ^k (r), where the latter reads in terms of the Pauli spinors 
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<p±(r) for both positronium configurations, cf. (13.191) : 

Wjo(r) = ( %(r-) = ifl{f) • <p + (r) + <pl{r) • <^_(f) . (4.1) 

Furthermore, if the decompositions (I3.51al) - (13.51bl) of the Pauli spinors with re- 
spect to the (double valued !) c<>Basis are inserted here, one finds the desired 
Maxwell density ^jo(r) appearing in terms of the (unique) wave amplitudes 
R±(r,$) and S±(r, i?) ( I3.53al) -( f3.53bl) in the following form 

(p) Jo(r) = R ?+f + S ++ S - . (4.2) 
47r r sm v 

However, it is not necessary to work with this exact form of the Maxwell charge 
density; but for our non-relativistic approximation we can make use of the fact 
that the spin-up and spin-down components of the time-independent Pauli spinors 
(I3.51al) -( f3.51bl) become decoupled (see the discussion of this point below equation 
(I3.57dp ). This enables us to restrict ourselves to the spin-up configurations alone 
(i.e. S± => 0), for which the charge density ^jo(r) (14. 2p reduces to 

6 2 I 6 2 

{p) m t±±i_ . ( 4.3) 

47T r sm v 

But here one may go one step further and restrict oneself to the non-relativistic 
limit which means putting i?_ to zero (see the discussion of this point below 
equation (13.61dl) ). Thus we ultimately arrive at the following approximative form 
of the one-particle charge density 

47T r sm v 

where the non-relativistic wave amplitude R+(r) R(r) has already been specified 
by equation (I3.66al) . The normalization of the non-relativistic density ^jo(f) (I4.4p 
is the following, in view of the fact that both matter fields Vv(r") (a = 1, 2) carry 
just one charge unit: 

1 = J d 3 f M j (r) = ^ J d 2 fR 2 (r) = ^ J drrR 2 (r) , (4.5) 
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with the two- and three-dimensional volume elements being given by 



d 3 f = r 2 sin i9 dr dfi 
d 2 r = r dr d$ . 



(4.6a) 
(4.6b) 



Here, the non-relativistic wave amplitude R(r) to be used hereafter for the calcu- 
lation of the positronium groundstate data is the following, cf. ( 13.66aj) 



R(r) = 



irr: 



exp 



r 
r* 



(4.7) 



The length parameter r* will subsequently be considered as a variational parameter 
for the minimalization of the groundstate energy. 



B. Spherically Symmetric Approximation 

But now that a suitable form of non-relativistic wave amplitude R(r) (14.71) is 
fixed, one can substitute this into the integral (13.581) for the determination of the 
electric potential ^Aq{t) which then appears in the following form: 

» M * = ?±jJ&wM^. (4.8) 

Observe that the asymptotic Coulomb form is guaranteed here by just the nor- 
malization condition (14. 5p . However, the present result (14.81) is not yet the desired 
form of the electric potential because it is not spherically symmetric. This obliges 
us to look for its spherically symmetric approximation which will be found by 
expanding the denominater of the integrand (14.81) as follows: 

1 1 f r*r' 3 / r • r ' \ 2 1 . 

1 • ~> ~r> ■ ^ +••• • (4-9) 



?"* — 1 1 "sj T^ 1 ~~\~ ^ ^ 2 T* 2 2 V V 2 ~\~ T* 2 

Clearly, such an expansion induces an analogous expansion of the electric potential 
H4 (r) g3D, i.e. one finds 

(p)A {r) = «A (r) + {II) A {r, 0) + . . . , (4.10) 
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with the spherically symmetric contribution ^Ao(r) being given by 
and correspondingly the first anisotropic correction ( /7 l4o(r, $) by 

WiM )^/ ffj/ MjM! . (4 , 2) 

o7T J ^ r 2 _|_ r /2j 2 

Naturally, one would like to neglect here the first anisotropic correction 
( 7/ l4o( r ? (and all higher-order terms) and retain only the spherically symmetric 
contribution "l4 (r) (14. lip . But this is allowed only if the anisotropic correction 
terms are actually much smaller than the spherically symmetric term. More con- 
cretely, it is easy to see that the first anisotropic correction (I4.12p is factorized 
according to 

( 7/ U (r, <&) = {II) A (r) • (// U (tf) (4.13) 
with the radial factor ^A (r) being given by 

^Uo(r) = a s r 2 [ dr' r' 3 • - 5 , (4.14) 



I _|_ ylA j 2 

and analogously the angular factor ^Ao^) by 

( 7/ U (tf) = J J dd'd^' (?' ■ r f = ^(1 + cos 2 0) . (4.15) 
(r = r/\\r\\ , etc.) 

Consequently, since this angular factor is of the order of unity, the anisotropic 
corrections could be safely neglected if the radial factor ^Ao(r) (14.141) is much 
smaller than the spherically symmetric potential ^Aq{t) (14.111) . It may be some- 
what difficult to show this for the most general situation, but a rough estimate 
can be carried out for a concrete model density R(r) 2 , i.e. the step function (see 
Appendix). Indeed for this special demonstration, the magnitude of the radial 
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factor ( 7/ l4o(r) amounts to (at most) a few percent of the spherically symmetric 
potential {I) A (r) fTCTTTl . see fie.l. 

Furthermore, the radial factor ( /7 l4 ( r ) (14.141) vanishes at the origin (r = 0) so 
that the exact potential ^A (r, $) (14. 8 p looks in any case spherically symmetric 
around the origin. And additionally it looks spherically symmetric at infinity 
(r — > oo) since it approaches the (spherically symmetric) Coulomb potential in that 
asymptotic region, cf. (13.581) . Thus the anisotropy of the electric potential ( p l4 (f) 
can be of some (minor) relevance only in the intermediate region (fig.l). In any 
case, the responsibility for the occurrence of binding is surely due to the spherically 
symmetric part ( 7 l4 (r) alone; and the effect of anisotropy will merely consist in 
some minor shift of the energy levels whose magnitude will be determined mainly 
by the spherically symmetric part ^Ao(r). For this reason, we will subsequently 
restrict ourselves to the use of that spherically symmetric approximation ^Aq(t) 
of the electric potential. 

C. Struve-Neumann Potential 

Clearly, even if one has decided to be satisfied with the spherically symmetric 
approximation ( 7 ^4 (r) (14. lip of the potential, one nevertheless wants to have this 
approximation as realistic as possible, i.e. one becomes faced now with the problem 
of an appropriate choice of the wave amplitude R(r) in equation (14. lip . Since in 
the asymptotic region (r — > oo) the approximative potential ^Ao(r) approaches 
the Coulomb potential, cf. (I3.58p . and since in this case the non-relativistic wave 
equation (I3.68P has the normalized wave amplitude R(r) (14.71) as its exact solution, 
one will now take this function (14.71) as a trial function and will calculate with its 
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help the associated trial potential ( J l4 (r) (14.111) 



d 2 r 



j iy> 2 

{Struve-N eumann potential) . 



as ( _2 \ 
7r Vr* 



(4.16) 



Indeed if one lets the length parameter r* tend to zero (r* — > 0), the trial charge 
distribution R(r) 2 shrinks to a point-like distribution whose associated potential 
^Ao(r) (I4.16P actually tends then to the Coulomb potential. However, the solution 
of the (non-relativistic) wave equation (13.681) is by no means point-like, but one 
rather expects that the realistic value of the ansatz parameter r* will fall in the 
vicinity of the Bohr radius a& ( 13. 61 aft ; and the tentative hypothesis is now that 
for this order of magnitude of r*, there will arise some potential t 1 ^o(r) from the 
prescription ( 14.161) so that its associated solution R(r) of the wave equation ( 13. 68 j) 
has an extension again of the order of magnitude of the starting value r* w ae- 
The question is now, how does such a potential [ 1 'y4 (r) look like? Especially, is 
there some similarity to the original Coulomb potential? 

In order to answer these questions, one simply calculates exactly the integral 
in ( I4.16P (see ref. 211 ]) and finds the following form: 



2x2 



dr 



r'=0 



f r' ■ exp [ — 2^-] 

j 2 



V r* / r 



(4.17) 



where the Struve-Neumann screening factor f$N is given by 



r* / V r* 



- 1 



(4.18) 



Here H\{z) is the first Struve function |21| 

(-!)»(§)* 



2 

7T 



7T 



EAz) 



(4.19) 



n=o ( n+ I).r(n + i)' 

and Ni (z) is the first Neumann function which admits a similar expansion as the 
Struve function. But the physically important point with the Struve-Neumann 
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potential (I4.17P is its asymptotic behaviour for r — > (origin) and r —>■ oo. For 
;he latter case, one has the following expansion for the first Weber function E\(z) 
211 



EAz) 



-NAz) - 



O 



1 



7XZ" \\Z\ 

Therefore the screening factor Jsn (14. 18ft becomes in the asymptotic region (\z\ 
oo) 



(4.20) 



fsN(z) 



7T 



- - E l {z)) - N x {z) 
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1 



E 1 (z)+N 1 (z) 



1+ °(w)- (4 - 21) 



But since the screening factor becomes thus unity far away from the origin, the 
Struve- Neumann potential ^^ (r) (I4.17P actually adopts the exact Coulomb form 
( I3.58P for r — ► oo. Clearly, this result is just what must be expected if the integral 
in (14.171) is correctly calculated; but the important point refers of course to the 
origin (r — > 0). 

For \z\ — > 0, the Struve and Neumann functions have the following behaviour 



H x (z) =o{\z\ 2 ) 

N 1 (z) = -- + -\n Z - + o(\z\) 

7TZ 71 2 



(4.22a) 
(4.22b) 



and therefore the Struve-Neumann screening factor fsN{ z ) (I4.18P becomes in the 
vicinity of the origin 



Jsn{z) 



z-z 2 - Z -\n Z - + o{\z\ z ) 



(4.23) 



Thus the screening factor weakens somewhat the Coulomb singularity at the origin, 
but does not eliminate it completely: The Struve-Neumann potential ' 1 ^ (r) (I4.17P 
behaves at the origin in the following way (r < r«) 



2as / 2r 1 /2r\ 2 , / r > , 
— (1 o - \n(-)+o(r 2 " 



(4.24) 
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This result says that the potential remains finite at the origin 

W4 (0) = ^ , (4.25) 

r* 

and this holds also for the radial component ^E r of the corresponding electric field 
W£ ( f) = -V [1 U (r), cf. (13351) : 

\im [1] E r ( y r) = a s (—) 2 , (4.26) 



see fig.l. 

Such a less singular behaviour at the origin (r = 0) has a very important 
consequence. The field energy being carried by the electrostatic field ^Aq( 
is given by the integral of the corresponding energy density as 

E$ = [ d 3 f llV^lyr)! 2 , (4.27) 



B3 



Aw as _ 

see also the energy functional below. However, if the present interaction po- 
tential ' 1 ^o(r) would exactly obey the Coulomb form with its singularity (^ 1/r) 
at the origin, the energy integral (I4.27P would diverge and thus the positronium 
groundstate energy would be infinite. Therefore, the partial regularization of the 
Coulomb interaction potential is necessary in order to equip the RST field config- 
uration with a well-defined energy content. But this point requires now a more 
thorough inspection. 



V. ENERGY FUNCTIONAL 



Quite generally speaking, the energy carried by a bound field configuration 
provides an ideal handle for testing the physical relevance of any theory because 
the atomic and molecular energy levels are immediately accessible to spectroscopic 
observation. For the present situation, the total energy Et of a bound RST field 
configuration appears as the sum of a matter part (Ed, say) and of a gauge field 
part (E G ), i.e. 

E T = E D + E G . (5.1) 
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All three kinds of energy are calculated as the spatial integrals of the corresponding 
energy densities T 00 (r), i.e. 



E T = J d 3 f ( T >T 00 (r) (5.2a) 
E D = J d 3 r^ D >T 00 (r) (5.2b) 
E G = J d 3 r^ G >T 00 (f) , (5.2c) 

where these energy densities are nothing else than the time components of the 
corresponding energy-momentum densities ^T^, ^^T^v and ^ G ^T^ U , resp. (see the 
extended discussions of this in the preceding papers). 



A. Gauge Field Energy Eq 

On the other hand, these energy-momentum densities T^ u are built up by the 
corresponding matter field (^) and gauge field (J 7 ^) so that, e.g., the electromag- 
netic gauge field energy Eq will appear as a sum of the energy content E^ due 
to the electric fields E a (r) (a = 1,2) and of the energy content E^ due to the 
magnetic fields H a (r): 

E G ^E R = 4 e) + 4 m) = J d 3 f {E^r) • E 2 (f) + H x {f) • H 2 (f)} , (5.3) 

cf. (14.271) . (For identical particles, the gauge field energy Eq contains also the 
exchange energy Eq, which however is missing here because positronium consists 
of two non-identical particles). 

An important point with the electric type of gauge field energy E$ Q refers 
now to the fact that this quantity can be recast also in an alternative form in 
terms of currents and potentials. To this end, one reconsiders the electrostatic 
Poisson equation (13.181) from which the following identity can be deduced, namely 
by multiplying through with the potential ^A (r) and integrating over, under use 
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of GauB' integral theorem: 

N G = J c/¥|rsin^||V (p U (r,^)|| 2 -as {p Uo(r)-(^ + /?_ 2 + S + 2 + ^ 2 )} =0. 

(5.4) 

Here, the specific relativistic form of the charge density ^7o(r) (14.21) has been used, 
which of course for our present non-relativistic purposes may be approximated as 
shown by equation ( 14.41) . As we shall readily see, this Poisson identity ( 15.41) turns 
out to be crucial for the subsequent principle of minimal energy. 



B. Matter Energy Ed 

The matter energy Ed is a somewhat more complicated object than the gauge 
field energy Eq. First, since the matter energy density ^ D ^Too{r) is built up by the 
time derivatives of the matter field cf. f 1 2 . 2 1 f) 

ihc 



*T (V V) - (V V)T * , (5.5) 



and since furthermore the time derivatives of \1/ (I3.1al) - (13.1bl) must explicitly con- 
tain the mass eigenvalues M a (=>- ±M#), the matter energy Ed (15.2bj) of the 
two-particle system must necessarily be built up by twice the one-particle mass 
eigenvalue M* 

E D = 2 M*c 2 - 2 M^c 2 , (5.6) 

where the mass equivalent (M^c 2 ) of the electrostatic gauge field energy ( I5.3P 
is defined through 



M fc 2 = ~jf d 2 f^A (r){Rl + R 2 + S 2 + S 2 } . 



(5.7) 



Indeed, the Poisson identity (I5.4f) ensures just the identity of the electrostatic gauge 
field energy E^ ( 15.31) and its mass equivalent M^c 2 (15.71) ! Now the reason, why 
the electrostatic interaction energy (or its mass equivalent, resp.) must be 
subtracted from the mass eigenvalue M*c 2 (for any particle) becomes immediately 
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clear: since the mass functional M*c 2 (see below) is composed of the rest mass 
energy (Mc 2 ), the kinetic energy (Tj^ n ) and the electrostatic interaction energy 
M^c 2 one has to subtract the latter (as a two-particle effect) from the mass 
functional M*c 2 in order to find the two-particle matter energy Ed 115. 6ft as the 
sum of pure one-particle contributions! Clearly it is very pleasant to see this 
subtraction mechanism working automatically by the right choice of the energy- 
momentum density ^■ D 'T (lv (12.211) for the matter field ^> . 

But with both constituents Ed and Eq of the total energy E? (15. ip being 
elaborated, one may now face the question whether perhaps this energy functional 
Et adopts its minimally possible value just for the solutions of the RST eigenvalue 
equations fl332iD - fl332dl) ? 

C. Action Principle 

The preceding RST dynamics, as a coupled set of matter and gauge field equa- 
tions, can also be deduced from a variational principle. The existence of such 
a principle is important for finding approximative (variational) solutions to the 
eigenvalue equations for those situations where it is hard to find the exact so- 
lutions. Indeed, we will subsequently present such a variational solution for the 
positronium groundstate. This will not only support the claim that exact solutions 
of the RST eigenvalue system do really exist, but it yields also valuable hints on 
further improvements of those variational techniques. 

Mostly, the action principles are based upon a Lagrangean density, here 
L R st , which is integrated over some space-time region in order to yield 
the action integral (Wrst) 



The extremalization of this action integral with respect to the matter field 
then reproduces the matter field equations (here the two-particle Dirac equation 




(5.8) 
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(I2.18P ); and similarly the variation of Wrst with respect to the bundle connection 
A^(x) yields the gauge field dynamics (i.e. the non-Abelian Maxwell equations 
(12.38!) ). For the present case of RST, the Lagrangean Lrst splits up into two 
parts, namely the matter part (Ld) and the gauge field part (Lq): 

AJ = £d[*] + Lo[A^} . (5.9) 

Here the matter part is given by [?| 

trie r — — 
£d[#] = — [^(V^) - (D^)r^J , (5.10) 

and the gauge field part looks as follows 

he 

LolA ' ] = m^ a K ^ F (5 - n) 

a, 6=1 

Observe here that the four- dimensional structure algebra u(2) is not sweeped out 
completely by the positronium curvature because the two-particle structure 
group U{2) becomes reduced to U'{2) = U(l) x U(l), see the discussion of the 
reduced Maxwell equations fl2.47ah - fl2.47bn . 

In order to be convinced that the present positronium eigenvalue systems 
( I3~52a1) -( i332dl) and (l3TM -( !3~56d1l actually do emerge as the Euler-Lagrange vari- 



ational equations from the claimed action integral (j5.8p - (j5.1ip . one substitutes the 
stationary ansatz (I3.1al) - (13.1bl) into the matter Lagrangean Ld (15.101) . uses the de- 
composition (13 .4p of the Dirac spinors into Pauli spinors together with the ansatz 
(I3.51al) -( l3.51bl) for those Pauli spinors, and thus finds the matter Lagrangean 
Ld^] being splitted up into three parts: kinetic term (L D %m ' ) ) and electric (L^) 
plus magnetic (L^) interaction term, i.e. 

L D m=L^ n) + L^ + L^ . (5.12) 
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Here the electromagnetic interaction terms are the simpler ones and are given by 

= 2 [M*c 2 + he {p) A {f)] ■ (%(r) (5.13a) 
L { ™ ] = ±2 he • {p/ %(r) , (5.13b) 

where the common mass eigenvalue M* is included in the electric part ( 15.13ajl : 
furthermore the upper /lower sign refers to the ortho/para case, and the common 
charge and current densities ^k (r) (I3.19P and ^^k^r) (I3.20p . (13.361) are given in 
terms of the wave amplitudes R±, S± by 

(%(f) = R l+f- + S l + S - (5 . 14a) 

{b/p) = sin^ [R+R_ - S+S-] - cos?? [S+R- + R+S-] 

^ ' 2irr sintf ' 1 ' } 

In contrast to these interaction terms of the matter part Ld (15.121) . its kinetic 
contribution looks somewhat more complicated but is the same for the ortho- 
and para-configurations: 



j {kin) _ he 



L> 2nr sin-# 



R (^^-l^L)-R ( dR + | ldS +\ 

+ dr r d$ dr r dd 

+ dr r d§ dr r d°& 



+ -(R + -R_ + S+- S\ - ^ (Rl - Rl + - S 2 )} . (5.15) 
r h 

But now that the matter Lagrangean fl 5 . 1 2 [) is explicitly known, one substi- 
tutes this into the action integral Wrst (15.81) . which itself splits up into two parts 
according to the splitting of the Lagrangean L RST (15. 9p 

W RST = W D + W G (5.16a) 
W D = J d 4 x L D [*] (5.16b) 

W G = J d A xL G [A n } , (5.16c) 
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and then the variational procedure for the matter action Wd f!5.16bj) with respect 
to the wave amplitudes R±, S± just yields the mass eigenvalue equations (I3.52al) - 
( I3.52dl) for the para-case and (I3.56al) -( f3.56dl) for the ortho-case. Here, there is a 
nice consistency check which, however, has a deeper meaning. Namely, the value of 
the matter action functional Wt> ( 15.16bl) upon the solutions of the mass eigenvalue 
equations must vanish because the latter are linear and homogeneous. Indeed, for 
any solution of such type of equations one can easily generate new solutions by 
simple muliplication by some constant c* (i.e. R± =>- c* ■ R±, S± =>- c* • S±). But 
this process of multiplication cannot change the value of the matter functional Wn 
( 15.16bl) since the multiplied wave amplitudes are also solutions of the variational 
procedure and do continously approach the original solution for the constant c* 
tending to unity (c* — > 1). And actually, substituting the derivatives of the wave 
amplitudes from the mass eigenvalue equations into the matter functional Wn 
( 15.16bl) lets this functional adopt the value of zero! This is the reason why the 
present action principle (5Wr> = 0) for the matter fields R±, S± does not require 
some constraint upon the latter fields (e.g. normalization condition)! 

This situation changes if one considers now the gauge fields, i.e. the electro- 
static potential ( p l4 (r) (its magnetic counterpart is not considered here because 
we restrict ourselves to the electrostatic approximation). The field equation for 
( p l4 (r) is the Poisson equation (13.181) which is to be deduced from the variational 
principle for the RST action integral (15.81) . Here it is sufficient to take into account 
only those parts W^ T of Hrst which contain the electric potential ( p l4 (r), i.e. 




RST — 



w b + W G 



(5.17) 



with 




(5.18b) 



(5.18a) 
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where the electric part of the matter Lagrangian Ld is given by equation 
( j5,13aj) ; and similarly the electric part of the gauge field Lagrangean L G ( 15.111) is 
given by 

4 e) = --^||V^ (f)|| 2 . (5.19) 

The important difference between the matter and gauge fields is now that the 
Poisson equation (13. 18ft is not homogeneous, and therefore one cannot multiply 
its solution WAo(r) by some constant c* in order to obtain a further solution 
c* • ^14.0 (f). It is just this peculiarity of the gauge potential ^A (r) which gen- 
erates the former Poisson constraint (15. 4p ! In order to become convinced of this, 
substitute the modified solution c*-^A (r) into the electric part W^ T ( 15.171) of the 
action integral so that this quantity becomes an ordinary function of the constant 
c*: VFftg T (c*). Since the original potential ^A (f) is presumed to extremalize the 
electric action H^g T (I5.17p . one is led to the following conclusion 



(e) 

0, (5.20) 



dc* 

and this is found to be nothing else than the Poisson identity ( 15. 4ft . Thus, through 
the present deduction, the Poisson identity may be interpreted to be a constraint 
for the variational procedure (see below for the use of this constraint in connection 
with the extremalization of the binding energy). 



D. Mass Functional 



Surely, it is a very pleasant feature of a theory if it can be based upon an action 
principle (here: 5Wrst = 0); but with reference to the treatment of bound systems 
the action integral W is not immediately linked to the binding energy which itself, 
on the other hand, is the quantity of physical interest. Therefore the question must 
necessarily arise whether, besides the action PVrst, the binding energy perhaps 
does also adopt an extremal value for the bound states (5Et = 0)? If so, how does 
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the corresponding energy principle look like and what is its relation to the action 
principle! The answer to these questions is most conveniently obtained by first 
considering the mass functional (Mtc 2 , say). 

For any solution of the mass eigenvalue equations (I3.52al) -( f3.52dl) or (I3.56al) - 
( 13.56dl) . resp., the mass eigenvalue M*c 2 adopts a well-defined real number, and this 
number can be represented in terms of the wave amplitudes R±, S± by multiplying 
through the eigenvalue equations by the associate amplitudes, integrating over 
whole three-space and adding up all four contributions. After the mass eigenvalue 
M*c 2 has been represented in this way, it may be reinterpreted also as the value 
of the corresponding functional M*c 2 =>- Mtc 2 [R±, S±] upon the solutions R±, S± 
of the mass eigenvalue problem. Thus the desired mass functional M^c 2 is found 



to appear in the following general form 

M T c 2 



Z 2 ■ Mc 2 + 2 (T r + T#) + M^c 2 + M { ™ ] c 2 . (5.21) 

Here, the meaning of the four constituents is the following: First, it must be 
mentioned that the constraint of wave function normalization has to be applied in 
its fully relativistic form 

1 = J d 3 f^j (r) = ^J d 2 r{R 2 + R 2 + S 2 + S 2 } , (5.22) 

whose non-relativistic approximation has already been specified by the preceding 
equation (14.51) . Next, the first term on the right-hand-side of equation fl 5 . 2 1 [) is a 
kind of modified rest mass energy (Mc 2 ) where the renormalization factor Z 2 is 
given by 



Z< 



d 2 r{R 2 -R 2 + S 2 -S 2 } . 



Furthermore, the relativistic form of the kinetic energy T^ in 
terms of the wave amplitudes R±, S± 



(5.23) 

T r + T$ reads in 



d 2 r 



he 
T 

^ / ■ - 



R. 



dR ± _R ± 8{rR. 
dr 
dS 4 



S 4 



dr 

d<JL 



+ R + 



dS + 
dr 
dS- 
~~dT 



S+ djrS.) 
r dr 
d(R + 



d$ 



(5.24a) 
(5.24b) 
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Finally, the last two terms are the mass-energy equivalents of the energy Er con- 
tained in the electromagnetic field modes, i.e. for the electric type 



Clearly, the upper/lower case of the latter equation refers again to the ortho/para- 
positronium. 

But with the mass functional MyC 2 being explicitly known, it is a standard 
variational procedure to actually deduce thereof the mass eigenvalue equations 
( j3.52a[) - (|3.52dp or (I3.56al) - fl3.56dl) . resp. Clearly, the relativistic normalization con- 
dition (15.221) has to be respected here, but this can be easily done via the method 
of Lagrangean multipliers (see below). The mass functional Mtc 2 (I5.2ip is not 
identical with the desired energy functional, but its closer inspection will help to 
construct the latter one. 

It seems very natural that the mass functional M^c 2 (15.211) consists of an electric 
and magnetic interaction term, and one expects that the other two terms should 
represent the rest mass and kinetic energies, resp. But why does the kinetic energy 
Tkin (= T r + T$) appear with a pre-factor of two, and why does there emerge a 
renormalization factor Z 2 in front of the rest mass term Mc 2 ? These questions 
are clarified most transparantly by passing to the non-relativistic limit of the mass 
functional. For this purpose, one eliminates the negative Pauli amplitudes £_ 
by expressing them approximately in terms of their positive counterparts R + , S + 




(5.25) 



and similarly for the magnetic type 
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via the mass eigenvalue equations, i.e. one puts 

R -~2Mc\ dr + r M J (5 ' 27a) 

§ - ~ 2iV^{^T " r^} ' (5 ' 27b) 
and if this is substituted into the relativistic kinetic energies T r (I5.24al) and T$ 
(15.24bl) one obtains their non-relativistic approximations ^E^m and ^E^m as 

T r =» <^,„ = £ / dV U^)' + ( d 4±) 2 ) (5.28a) 



4M J I V dr J V dr 
T^m EHn = ^ [d*rta a 4±)\a a M 2 } (5.28b) 



4M J IVr d'd J \r d'd 

h 2 ff , , n r&R + 95+ (AS',_ f 7/? 4 
+ ttt / / dr dv 



AM J J " ~" L dr dtf <9r 
Thus the non-relativistic approximation i?/^ (= ^E^ + ^E^) of the relativistic 
kinetic energy T kin (= T r + T&) is found to appear in the following form 

Ekm = W j d2 r{\\VR+\\ 2 +\\VS + f}+E w . (5.29) 

Here the "winding energy" E\y is due to the map from the (r, i?)-plane into the 
(R+, S+)-space and is defined (as usual for winding numbers) in terms of the 
corresponding Jacobian as 

E * * m II " rM • < 5 - 3 °> 

Indeed, apart from the emergence of this winding energy (to be neglected for the 

present purposes), the kinetic energy E kin (15.291) is just the result to be expected 

from the non-relativistic approximation of the mass eigenvalue equations 

fi 2 / d 2 1 d 1 d 2 \ ~ 

+ -^- + ^-^)R-hc^A -R = E s -R (5.31a) 



2M \dr 2 r dr r 2 d'd 2 

-If + \l + hw)~ s ~ ^' ,A ° - § - Es - § (5 ' 31b) 

Obviously, the spin-up (~ R) and spin-down (~ S) components of the non- 
relativistic approximation are decoupled, because the magnetic interactions are 
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neglected (see equation (13.681) for the spherically symmetric approximation hereof). 
Concentrating here upon the spin-up solutions (15.31a|) . one finds the Schrodinger 
energy functional Es[R] by simply multiplying through that equation by R and 
integrating by parts as 

E ^mj d f { ( * ) + 5* ( m) } - T / ** P>A ° ■ R • (5 32) 

where the non-relativistic normalization condition has been applied, cf. (14. 5ft 

d 2 r R{r,$) 2 = 2 . (5.33) 



But on the other hand, this Schrodinger energy functional Eg (15.321) must turn 
out as the non-relativistic electrostatic approximation of the relativistic mass func- 
tional Mj-c 2 (I5.2ip . where the non-relativistic form of the electric mass-energy 
equivalent M^c 2 (15.251) in equation (I5.32p is self-evident. However, the relativis- 
tic kinetic energy T^ in (= T r + T$) appears with a factor of two in the relativistic 
form (I5.2ip and, on the other hand, appears in the usual way in the non-relativistic 
form (I5.32j) . cf. (15.291) . In order to clarify this final point, one subtracts the rest- 
mass energy Mc 2 from the relativistic mass functional Mtc 2 in order to obtain the 
non-relativistic Schrodinger functional E$ 

E s = M T c 2 - Mc 2 , (5.34) 

and furthermore one looks also for the non-relativistic approximation of the renor- 
malization factor Z 2 (I5.23P by use of the approximations (l5.27aD - fl5.27bl) which 
yields 

Z " = 1 - » ( 5 - 35 ) 
with the non-relativistic E^in being given by equation (I5.29p apart from the winding 
energy Ew- Consequently if this is substituted back into the relativistic mass 
functional M T c 2 flCTj) . one finally gets by observation of the non-relativistic form 
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Ekin (I5.29P and with the neglection of the magnetic contribution M^ 1 c 2 (15.261) for 
the Schrodinger energy functional Eg (15.341) : 



Thus the non-relativistic approximation (I5.36P of the relativistic mass functional 
Mt(? yields exactly the Schrodinger energy functional Eg (15. 32ft ; and this of course 
supports the confidence into the mass functional M^c 2 as the correct relativistic 
generalization of the non-relativistic Schrodinger functional Eg. 

This viewpoint is even further supported by considering the variational equa- 
tions due to both energy/mass functionals. Namely, it is well-known that the 
ordinary non-relativistic Schrodinger equation for a fixed potential ^A arises by 
extremalization of the functional Eg under the constraint that the wave func- 



any textbook about quantum mechanics, e.g. ref.s[16(,|l7|] for the Ritz variational 



method. But if the present RST eigenvalue systems (I3.52a[) - (13.52dj) . or (I3.56al) - 
(13.56dl) . resp., do really represent the correct relativistic generalization of the non- 
relativistic Schrodinger equation, then the general logic demands that these rela- 
tivistic mass eigenvalue equations must appear as the variational equations due to 
the mass functional Mtc 2 in the same sense as the ordinary Schrodinger equation 
arises as the variational (i.e. Euler-Lagrange) equation due to the non-relativistic 
Schrodinger functional Eg\ And indeed, if the corresponding variational proce- 
dure for the relativistic functional M^c 2 (I5.2ip is carried through with regard of 
the original relativistic normalization condition (I5.22p . one arrives just at those 
RST eigenvalue systems for the ortho- and para-cases, resp. To be more specific, 
multiply the normalization condition (15.221) 



Eg = {Z 2 - l)Mc 2 + 2 E km + Mfc 2 = E kin + Mf c 2 



(5.36) 





present situation), see 




(5.37) 
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by some Lagrangean multiplier and add this to the mass functional M^c 2 (I5.2ip 
in order to obtain the modified functional M^c 2 

M T c 2 = M T c 2 + \ D ■ N D . (5.38) 

Now carry through the variational procedure for the modified functional MtC 2 
and find the RST eigenvalue systems such as f)3.52ap - fl3.52dj) with merely the mass 



eigenvalue M* being replaced by the Lagrangean multiplier A^, i.e. 

M,c 2 => -X D . (5.39) 

Summarizing, one can be sure now that one has succeeded in constructing the 
correct relativistic generalization of the ordinary, non-relativistic Schrodinger ap- 
proach; but one further generalizing step has to be done, and this refers to the 
interaction potential between the two particles: wheras the ordinary Schrodinger 
approach is based upon the instantaneous Coulomb potential (see the Schrodinger 
Hamiltonian H (12.221) ). the relativistic RST interaction potential ^Aq is a dynam- 
ical object which has to obey its own field equation (i.e. the Maxwell equation 
(I2.38p in the general case)! In this sense, RST is a more fundamental framework 
than the ordinary Schrodinger quantum mechanics. But more concretely, this 
means now that we finally have to generalize the preliminary functional Mj>c 2 
(I5.38P in such a way that also the Poisson equation (I3.18P for the electric poten- 
tial ^A Q arises as an Euler-Lagrange equation due to that desired more general 
functional. 



E. Principle of Minimal Energy 

Surely, it will not come as a surprise that the functional, which we are after, 
must have something to do with the total energy Et (15.11) . Therefore it is instruc- 
tive to inspect this functional in greater detail, where the magnetic interactions 
may be neglected again for the present purposes. 
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First, observe here that the electrostatic gauge field energy can be 

expressed exclusively in terms of the gradient of the electric potential ^A (r) 

m = [ d 3 f\\V (p) A (f)f , (5.40) 

47ra s J 

because the electric field strength E p appears as a gradient field whenever the 
exchange interactions are missing, i.e. (cf. (I3.15P ) 

E p (r) = -V (p U (f) (5.41) 

(for the generalization of this relationship in the presence of exchange interactions 
see ref. 22J). Next, consider the contribution of the matter energy Ed (15. 6p which 
also contains the electric potential ( p l4 (r), but in a form different from the mass 
equivalent M^c 2 ( 15. 7ft . Furthermore, the electric potential enters also the mass 
functional Mtc 2 (15.211) again in form of the mass equivalent M^c 2 . Therefore, 
in order to make the appearance of the electric potential in the energy functional 
Et (15.11) explicit, one writes down the latter functional (under use of the mass 
functional M T c 2 (ET2Ijl ) as 

E T = E D + E G = 2(M T c 2 - Mfc 2 ) + E^ ] + 

= 2Z 2 -Mc 2 + 4 T km + E^ + E^ + 2 M^c 2 . (5.42) 

But this, indeed, is a rather amazing result because Et contains now the electric 
potential ^Ao(r) exclusively in form of the gauge field energy E^ (E3QD, whereas 
the mass functional M-^c 2 ( I5.2ip contains the potential < - p ^ (r) exclusively in form 
of the mass equivalent c 2 ( 15.251) ! And this implies that neither the mass func- 
tional Mt ( 15.211) nor the energy functional Et (15.421) would lead us to the electric 
Poisson equation ( 13.181) as the Euler-Lagrange equation of the corresponding vari- 
ational procedure. 

Evidently, in order to escape from this dilemma, one has to put in here some new 
idea. This is done now by imposing a second constraint, besides the normalization 
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condition (I5.37p . namely the Poisson identity (15 .4p . Consequently we multiply 
both constraints (15.41) and (I5.37P with the Lagrangean multipliers X G and A^, 
resp., and add this to the original energy functional Et (15.421) in order to thus 
arrive at the modified energy functional Et'- 

Et = Et + 2 Xd ■ Ntj + Xg • Ng 

= 2Z 2 - Mc 2 + AT kin + E$ + 2 X D ■ N D + X G ■ N G . (5.43) 

Here the renormalization factor Z 2 is given by (15.231) . the kinetic energy T kin 
(= T r + T#) by (Q4^ - (jK24bjl . and finally the gauge field energy E^ by (l5~40l) . 
Moreover, the magnetic terms are omitted because we are satisfied for the mo- 
ment with the electrostatic approximation. However, carrying through now the 
variational procedure for the final result Et (15.431) actually yields just the (fully 
relativistic) eigenvalue equations (I3.52a|) - (13.52d|) or (l3.56aD - (13.56dl) . resp., with the 



first Lagrangean multiplier being given as before by equation (I5.39P and the 
second one (Xc) by 

X G = — ■ (5.44) 

a s 

And additionally, the variational procedure for E T (I5.43P with respect to the elec- 
tric potential ^A (f) (and magnetic potential ^14^ (r)) yields the desired Poisson 
equations, i.e. the electric one (13.181) and also its magnetic counterpart (to be sup- 
pressed here). 

Summarizing, it is possible to convert the general Hamiltonian-Lagrange action 
principle (^W-rst = 0) t° an energy minimalization principle (5Et = 0) for the 
bound states. The corresponding Euler-Lagrange equations are just the relativistic 
mass eigenvalue equations plus the Poisson equations for the interaction potentials 
of electric and magnetic type. Clearly, such a principle of minimal energy for 
bound states provides now a powerful approximation technique when looking for 
the binding energy. We will readily exploit this pleasant result, but will be also 
satisfied with the corresponding non-relativistic approximations. To this end it 
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is merely necassary to resort to the non-relativistic approximation of those terms 
emerging in the relativistic energy functional E? (15.431) . 

First, the non-relativistic forms of the renormalization factor Z 2 and the ki- 
netic energy T have already been specified by equations (15. 29ft and ( 15.351) . resp., 
where the spin-down component S+ is put to zero so that the winding energy Eyy 
(I5.30p does also vanish. Second, the constraint term Np due to the normalization 
condition is omitted completely (also in its non-relativistic form (I5.33P ) because as 
a possible trial function R(r) we resort to the type (14.71) which is a priori normal- 
ized to unity. However, in contrast to this, the Poisson constraint term Nq must 
be retained in any case. The reason is that if one takes the value of the wanted 
functional E$ upon some trial field configuration, which is not an exact solution 
of the corresponding variational equations, then the Poisson constraint keeps the 
(naturally inaccurate) trial value of E$ in the neighborhood of its proper mini- 
mum (see the example below). Thus the non-relativistic version E$ of the energy 
functional Et (I5.43[) is ultimately found as 

E s = 2 E km + E%> + 2 (A4 e) c 2 - 4 £) ) (5-45) 

Here it is very instructive to consider the difference between the conventional 
Schrodinger one-particle functional E$ ( 15.361) and the present two-particle RST 
result E$ (15.451) . It is true, the kinetic two-particle energy {2E kin ) is twice the 
conventional one-particle counterpart as displayed in equation (I5.36p . however the 
electrostatic interaction energy is described in a radically different way: whereas 
the conventional Schrodinger approach Es ( I5.36P relies here exclusively upon the 
mass equivalent M^c 2 (I5.25P in its non-relativistic form, the RST counterpart Es 
( I5.45P describes the interaction energy in terms of both the field energy E^r and its 
mass equivalent M^c 2 \ Clearly, for an exact solution of the RST energy eigenvalue 
problem both contributions E^ and M^c 2 are identical (on account of the Poisson 
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identity) so that they cancel for the non-relativistic RST functional Eg (I5.45j) 
which then adopts the conventional Hartree-Fock form. However for (non-exact) 
trial configurations, the difference of the field energy 4 e) and its mass equivalent 
M^c 2 becomes important for the degree of accuracy of the corresponding trial 
value of Es- This is the origin of the fact that RST predictions have the potential 
to come closer to the experimental numbers than the conventional Hartree-Fock 
approach. Clearly, such a pleasant result must be exemplified now by a brief 
numerical treatment of the positronium groundstate. 



VI. NUMERICAL DEMONSTRATION 



Being satisfied for the moment with the non-relativistic approximation, one 
inserts the trial function R(r) (14. 7p together with the Struve-Neumann potential 
t 1 U (r) fl4.17l) - fl4.18p into the non-relativistic form Es ( I5.43P of the energy func- 
tional which then becomes an ordinary function of the length paramter r* (^(r*), 
say). The positronium ground-state energy Eq emerges then as the minimal value 
of the function Es(r*) which may occur for r* = ro- 



Eq = E s 
dE s {r*) 



(6.1a) 

. (6.1b) 



This RST groundstate energy E must then be compared to the corresponding 
groundstate energy E HS predicted by the conventional Hartree-Schrodinger ap- 
proach in order to test which one of both approaches yields the more accurate 
numbers in comparison with the observational data. 



A. RST Prediction 



Concerning the non-relativistic RST functional Es (15.451) . it is highly informa- 
tive to consider any of its contributions separately. First, the one-particle kinetic 
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energy E^ n is taken from (I5.29P or (15.361) and is always of the general form 

E kin = ^— I ' d 2 r (^^I) = 2 sun al Mc 2 — 
AM J \ dr J s y, 

(y* = r*/a B ) • (6.2) 

Here e^m is a real number which is associated with the exponential trial function 
R(r) (S3]), through 

oo 

ekin = J dyy e~ 2y = ^ 
o 

(y = r/r„) , (6.3) 

and the typical atomic energy ("atomic unit", a.u.) is given by 

e 2 

a\ Mc 2 = — = 27,21 [eV] . (6.4) 

CLb 

Thus the one-particle energy E kin (16. 2p appears as 

13,61 



J kin 



y* 



eV] . (6.5) 



Next, the electrostatic field energy E ( k ] CTJD is found in a similar way by use 
of the Struve-Neumann potential [1 U (r) ([CT7] - (j4TTgl) as 



Eg = \ d 6 r , uv/ = -2 i% a 2 s Mc 2 ■ - . (6.6) 

The real number is given in terms of the dimensionless Struve-Neumann screen- 
ing factor fsN( x ) (I4.18P as 



00 n 2 



6^ — / dec X 



dx\ x 



0,3822 . (6.7) 



o 



Thus the field energy is found as 



-(e) 20,80 

~ [eV\ . (6 J 



69 



And finally, the mass equivalent M^c 2 (15.251) is found to appear in its non- 
relativistic form as follows: 

M fc 2 ~ f d 2 f ^A (r)R(r) 2 = -2^ a 2 s Mc 2 ■ - , (6.9) 

^ J y* 

with the real number /i^ being given by 

oo 

$ = j dx e~ x f SN (x) ~ 0,4348 , (6.10) 
o 

which then yields for the desired mass equivalent 

M (e) c 2 „ [eV] . (6.11) 

y* 

Observe here that the kinetic energy Eu n (~ y* 2 ) is positive, whereas the 
electric contributions (~ y~ l ) and M^c 2 (~ y^ 1 ) are negative so that the 
expected equilibrium point (16.1bl) of the energy E s {y*) fl5.43j) 
~ 27,21 20,80 / 23,66 20,80x 

Es[y*) = — 2 + 2 ' + ) eV 4 

y* v y* y* 1 



27,21 26,52 

y* y* 



[eV] (6.12) 



does really exist and is found as 



yo==^= ^ R * 1,95 , (6.13) 



see fig. 2. Therefore the corresponding groundstate energy E Q (16. lap of positronium 
is obtained as 

E = -6,48 eV . (6.14) 

This differs from the experimental value (~ 6,80 [eV]) by 0,32 [eV] and thus falls 
within the error estimate (~ 1 [eV]) of the Appendix. The important point here 
is that the bracket term in the energy functionals (15.451) and (16.121) is an essential 
element of our present RST treatment of the positronium groundstate because it 
contributes roughly 3 [eV] (~ 50%) to the binding energy. Since however this term 
is missing in the conventional Hartree-Schrodinger approach, one expects that the 
latter approach will predict the groundstate energy with an error of (roughly) 50%! 

70 



B. Hartree-Schrodinger Prediction 



According to the general belief, the conventional Schrodinger equation 

r 2 ) = {S) E V(f u f 2 ) (6.15) 

would yield the exact groundstate energy ^E (apart from relativistic and QED 
effects) if this eigenvalue problem could be solved exactly. However, these exact 
solutions are not known in general (apart from such special cases as (12.221) ); and 
therefore one has to resort to certain approximation techniques, e.g. the Hartree- 
Fock method 16j, ll7j. But since the spin effects are neglected for the present 
purposes, that method reduces to the Hartree approach which is formally based 
upon the conventional one-particle Schrodinger equation ("Hartree-Schrodinger 
approach"). The essential point with this approximation method refers to try- 
ing a simple product ansatz for the two-particle wave function ^(ft,f 2 ) (16. 15ft in 
terms of one-particle wave functions ip(ft) and ip(f 2 ) 1 i-e. one tries for the present 
groundstate problem 

*(ri, f 2 ) = if} (ft) ■ ip(f 2 ) . (6.16) 



This ansatz may be inserted into the conventional action integral W$ 16l. Il7] 



W s = JJ d 3 ftd 3 f 2 ^*(ft,f 2 ) H*(r u r 2 ) , (6.17) 

which adopts its stationary values (5Ws = 0), within the set of normalized wave 
functions 

NS=¥ II ^ ' ^ fl ' f2 ) ~ 1 = ' (° 8 ) 

just for the solutions of the Schrodinger eigenvalue problem ( 16.151) . However for 
the trial functions (I6.16p . the Schrodinger action (I6.17P becomes converted to its 
Hartree-Schrodinger form Whs which looks as follows: 

W H s = 2 {s) E km + V H s , (6.19) 
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where the one-particle kinetic energy ^Ep-in is given as usual by 

{S) E km = ^J <P? W>*(r) • VV(r) (6.20) 
and the electrostatic interaction energy Vhs by 



v as = w*>M*ft)i' , (6 . 21) 

Ifi — r 2 | 

Obviously, the Hartree-Schrodinger action Whs (16.191) represents the energy of 
the two-particle system and therefore is the conventional counterpart of our non- 
relativistic RST energy functional Et (I5.43H . Therefore the corresponding energy 
minimalization due to the method of Lagrangean multipliers {S(Whs+^s'^s) — 0) 
yields the following one-particle eigenvalue equation of the Schrodinger type: 

- ^ A V^) - hc {s %(r) ■ i;(r) = -A s • ifrif) . (6.22) 

Here the Lagrangean multiplier As plays the role of the energy eigenvalue and the 
electric interaction potential ^A (r) is given in terms of the wave function ip(r) 
through 

(^o(f) = aJ d*r> . (6.23) 

It is true, both the energy eigenvalue equation f!6.22|) and the interaction potential 
(I6.23P look very similar to their non-relativistic RST counterparts, see equations 
(I5.31al) - (15.31bl) and (13.581) . resp. But a closer inspection will readily reveal very 
important differences. 

First, observe that the conventional eigenvalue equation (16.221) is based upon 
the ordinary Laplacean A in spherical polar coordinates (see (13.691) for its radial 
part), whereas the corresponding RST eigenvalue equation (I3.68P relies upon the 
spherical polar coordinate r as a cylindrical variable. Correspondingly, the RST 
normalization condition (14.51) for the wave function R(r) uses a two-dimensional 
integral, whereas the conventional formalism refers to a three-dimensional integra- 
tion, i.e. 

d 6 r \ip{r)V = 1 . (6.24) 
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For instance, if we adopt again an exponential trial function with length parameter 

r s 

^r) = -L=e- r ^ , (6.25) 
V 7rr £ 

which obeys the conventional normalization (I6.24p . this trial function must nec- 
essarily differ from its RST counerpart R(r) (J4TTJ) ! But nevertheless, the kinetic 
energies turn out to be identical, i.e. for the conventional case ( 16.201) one finds 

z iJs Us 

(y s ±r s /a B ) , (6.26) 



which actually is the same as its RST counterpart fl6.2p - fl6.3p . 

However, the crucial difference between both approaches does refer to the elec- 
trostatic interaction potentials (5) A)(f) fl6~2"3"]) vs. [1 U (0 f]CT7j) - (jl~Tg]) . For the 
conventional case 06.231) one finds by use of the trial wave function (16.2 



[^U (r) = ( 1 - e- 2r ^) ^ - ^ e" 2 ^ , (6.27) 
V / r r s 

and this is of a different type in comparison to the RST case ^A (r), see fig.l: 
the conventional potentials of the form (I6.23P have vanishing field strength ^E r 
(= d^A^iy) / dr) at the origin (r = 0), whereas the corresponding RST field 
strengths adopt a non- vanishing value, cf. (I4.26p ! Of course, this circumstance 
must have its consequences for the interaction energy Vhs (16.211) concentrated in 
the electrostatic field: 

V HS = -he j d 3 f ( 5 U (r)^(r) 2 = -^ajMc 2 ■ — ~ _ill2£ [ e y], (6.28) 

J o Us Us 

and this is now essentially different from the corresponding RST interaction energy 

dS3SD 

V RST = 2 i\4 e) c 2 - 4 e) ~ [eV] , (6.29) 

y* 
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cf. (I6.12p . Thus the conventional Hartree-Schrodinger energy Whs fj£>. 19|) is finally 
found as 

27 21 17f)f) 

W HS {y s ) = ^# - ^ [eV] , (6.30) 

Vs Vs 

in contrast to the RST case Es(y*) (I6.12p . Whereas in the latter case the RST 
energy adopts its minimal value ( 16.14ft for yo = 1,95, the minimal value of the con- 
ventional energy Whs (I6.30P occurs for y = 3,2 and is found as E ~ —2,65 [eV], 
see fig. 2. Clearly this is a very poor result when compared to the corresponding 
RST prediction of 6,48 [eV] (16. 14}) and to the experimental value of 6,80 [eV], resp., 
see fig. 2. 

Summarizing on the basis of the present groundstate calculations, it appears 
natural to assume that in the general situation RST will be superior to the Hartree- 
Fock approach, as far as the numerical coincidence of the theoretical predictions 
and the experimental values is concerned. Thus as the next task, there remains to 
be settled the (non-relativistic) competition between the conventional Schrodinger 
approach and RST. Here, as the Appendix demonstrates, the anisotropy of the 
RST field configurations must be taken into account. Indeed, such a competition 
seems to be of profound philosophical relevance, for if RST could predict the same 
numbers as the standard Schrodinger theory, the latter approach would lose its 
status of uniqueness! 



APPENDIX: SPHERICALLY SYMMETRIC APPROXIMATION 

In order to get some feeling of the extent to which the first anisotropic correction 
( 7/ l4o(r, $) (14.121) is smaller than the spherically symmetric contribution ^Ao(r) 
(14. lip , it is very instructive to consider a typical example, which admits to exactly 
calculate all the terms in question. 

Here, as the prototype of a well-localized charge distribution R(r) 2 , one may 
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adopt the following step function 

( 4 

< r < r. 



R\r) => ^R 2 (r) = { n ' r c ~ (A.l) 

, r > r c . 

It should be obvious that this (non-relativistic) charge density actually obeys the 
normalization condition (14.51) for any fixed length parameter r c . Furthermore, 
substituting this step density into the spherically symmetric approximation WAo(r) 
(14. lip yields the following model potential t c i4o(r): 

(/ U (r) l c k (r) = ^{y^T^-r} . (A.2) 

' c 

Of course, this potential ' c l4o(r) adopts the required Coulomb form ( 13. 58ft in 
the asymptotic region (r ^> r c ), see fig.l. If the length parameter r c is chosen to 
be infinitesimally small (r c — > 0), the Coulomb form does fill the whole three-space 
(0 < r < oo). The screening factor (f c , say) is found for general value of r c as 



/c (r)=2^Ul+(^) 2 -^} (A.3) 

qand thus is seen to tend to unity (/ c (oo) = 1) in the asymptotic region (r 3> r c ). 
On the other hand, the screening factor / c tends to zero at the origin (r — > 0); and 
therefore the model potential ' c lAo(r) ( 1A.2I) remains finite for r — >• 0, i.e. 

Uo(0) = ^ , (A.4) 



[c 



cf. the analogous behaviour of the Struve-Neumann potential MAj^O) (I4.25p . Ob- 
viously these are the common features of all those interaction potentials ^Ao(r) 
(14. lip due to the exotic states! (see fig.l). 

Next one wishes to check for the chosen step density (lA.ip . to what extent 
the first anisotropy correction ^A (r) (14.141) really is smaller than the spherically 
symmetric potential t c l4 (r) (IA.2[) . For this purpose, one substitutes the step 
density into the radial factor (I4.14p and thus finds 

mMr) ^ Mr) = ^{ r -^.^±^). (A . 5) 
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And indeed, this anisotropic correction vanishes at the origin (r = 0); and at 
infinity (r — > oo) it tends to zero faster (~ r -3 ) than the Coulomb potential 
(~ r -1 ), see fig.l. Therefore the neglection of the first anisotropic correction 
seems to be justified. 

However, perhaps a better measure of the relative magnitudes consists in the 
field energy (15.401) . Substituting there the decomposition (14.101) of the elec- 
trostatic potential ^A (r) into the radial part < - 7 ^ (r) and in its anisotropic cor- 
rections, one finds a similar decomposition of the electrostatic energy E^\ i.e 

E$ = + W£M + . . . (A.6) 

where the first energy contribution is due to the spherically-symmetric po- 

tential (/ l4 (r) 

W£fW = _^£_ / d 3 r ||V«A (r)f , (A.7) 
47ra s J 

and similarly the second contribution is built up by the spherically-symmetric 
potential together with the first anisotropic correction 

WgW = _ J d z r V (/ U (r) • V (// U (f) . (A.8) 

But the point here is now that both potentials ^A (r) and ^Ao(r) can be deter- 
mined exactly for the step function density, see (1A.2j) and (1A.5j) ; and therefore the 
corresponding energy integrals (A.7) and (A.8) are also known exactly, namely 

«4 e) [c] 4 e) = -(2tt - y)^ ~ -0,9499 t- (A.9a) 

= _M _ 137T £ 
fl fl 3 V 3 32 ; r c r c V 7 

Thus the first anisotropic energy correction (A. 9b) amounts to (roughly) 16% 
of the spherically symmetric energy (A. 9a)! This says that the restriction to the 
spherically-symmetric approximation lets us expect an uncertainty of our positro- 
nium groundstate result up to leV near the experimental value (of 6.8 eV). 
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Fig. 1: Spherically Symmetric Approximation 

The Struve-Neumann potential ^A (r) (4.17) (solid line) has a finite value 

(4.25) at the origin (r = 0), as well as the corresponding field strength ^E r (r) 

(4.26) . The Coulomb form ~ 1/r (3.58) is adopted at spatial infinity (r ~» oo). 
These features are common to all potentials due to a charge distribution of the 
type (4.2), e.g. also for the potential model ^A (r) (A. 2) generated by the cut-off 
charge (A.l) (broken line). The first anisotropic correction ^A (r) (A. 5) (in- 
tersected line) is considerably smaller than the spherically symmetric approxi- 
mation ^A (r) (A. 2) and is therefore neglected for the present rough estimate. 
The conventional Hartree-Schrodinger potential ^Ao(r) (6.27) (dotted line) has 
vanishing field strength at the origin, in contrast to the RST potentials ' 1 U (r) 
and l c U (r). 77 




Fig. 2: RST Energy Eg (6.12) and Hartree-Schrodinger 
Energy W m (6.30) 

The groundstate energy E (6.1a) is obtained from the minimum (6.1b) of the 
energy curves. RST and the conventional Hartree-Schrodinger approach agree in 
the kinetic energies but differ in the interaction energies Vhs and Vrst, see 
equations (6.28)-(6.29). This is the reason why the RST prediction (-6,48 [eV]) for 
the groundstate energy is closer to the experimental value (-6.80 [eV]) than the 
Hartree-Schrodinger prediction (-2.65 [eV]). 
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